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MODES  OF  CURRENT  CONDUCTION 

The  electrical  characteristics  of  solids  are  generally  classified  as 
representative  of  a  metal,  semiconductor,  or  insulator. 

A  metal  conducts  current  in  a  linear  manner  characterized  by  Ohm’s  Law. 
Semiconductors  and  insulators  also  follow  a  linear  mode  of  current  conduction 
if  the  applied  electric  field  is  small  and  the  external  and  internal  temperature 
is  moderate  (around  room  temperature).  At  higher  electric  fields  or  extremes 
of  temperature,  semiconductors  and  insulators  exhibit  nonlinear  electrical 
characteristics  representative  of. the  various  scattering,  generation,  and 
recombination  mechanisms  involved  in  the  current  conduction  process.  Oftentimes 
these  nonlinear  electrical  characteristics  can  give  an  insulator  the  low  field, 
moderate  temperature  characteristics  of  a  semiconductor.  Induced  nonlinear 
electrical  characteristics  can  also  give  a  semiconductor  the  low  field,  moder¬ 
ate  temperature  aspects  of  a  metal  or  an  insulator. 

Examples  of  this  type  of  nonlinear  behavior  are  quite  common  in  literature 
(Refs.  1,"2,  and  3).  One  common  example  is  the  ability  of  lightly  doped  or 
pure  semiconductors  to  resemble  insulators  electrically.  At  very  low  tempera¬ 
ture  at  the  opposite  extreme,  semiconductors  begin  to  possess  metallic-like 
values  of  conductivity  (Ref.  1).  High  temperatures  also  affect  some  insulators 
by  increasing  their  conductivity  to  an  extent.  Sometimes,  this  increase  in 
conductivity  is  enough  to  make  the  insulator  resemble  a  lightly  doped  semicon¬ 
ductor-electrical  ly. 

The  impression  of  high  electric  fields  (AC  or  DC)  across  semiconductors 
and  insulators  can  also  induce  changes  in  these  solids'  electrical  characteris¬ 
tics.  For  example,  high  fields  impressed  on  some  insulators  can  cause  permanent 
or  transitory  change  to  the  materials  electrical  characteristics.  Although  this 
is  accomplished  through  such  mechanisms  as  breakdown  in  quartz,  other  insulators 
(such  as  silicon  carbide)  resemble  a  semiconductor  at  high  fields  near  breakdown 
(Ref.  4).  Under  high  fields,  semiconductors  can  exhibit  metallic-like  behavior. 
In  particular,  this  type  of  behavior  is  exemplified  by  the  phenomenon  called 
second  breakdown  (Ref.  5).  The  phenomenon  known  as  second  breakdown  is  an 
unstable  transition  from  single-carrier,  space-charge-limited  current  to 
double-carrier  injection. 
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Alteration  of  the  electrical  characteristics  of  solids  can  also  be  attained 
by  chemical  means.  Doping  semiconductors  and  insulators  can  make  the  electrical 
characteristics  appear  more  metallic.  Generally  speaking,  chemically  altering 
the  electrical  characteristics  of  materials  is  somewhat  limited  unless  augmented 
by  temperature  or  electrical  changes. 

This  aspect  of  solids  (i.e.,  the  ability  to  change  their  electrical  char- 
I  acteristics  under  temperature  extremes  or  changes,  high  electric  fields,  or 

chemical  means)  demonstrates  that  the  nonlinear  aspects  of  solids  can  be  approx¬ 
imated  by  a  study  of  their  transition  from  one  class  of  solids  to  another.  In 
other  words,  the  nonlinear  behavior  of  say,  a  semiconductor,  can  be  modeled  as 
transitions  between  semiconductor-like  states  to  metallic-  or  insulator-1  ike 
states.  This  interpretation  of  nonlinear  phenomenon  in  solids  forms  the 
phenomenological  basis  for  the  analysis  that  follows.  From  it,  the  importance 
of  space-charge-limited  single-carrier  currents  and  double-carrier  currents 
will  become  obvious  in  explaining  observed  nonlinear  phenomena.  In  particular, 

I  this  provides  the  basis  for  explaining  the  preconditions  under  which  electrical¬ 

ly  unstable  states  of  operation  occur,  thus  resulting  in  device  transition  to 
second  breakdown.  « 
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SINGLE-CARRIER  CURRENT  FLOW 


The  electrical  characteristics  of  an  ideal  diode  show  a  high  degree  of 
conductivity  under  forward  .bias  with  a  small  voltage  drop  across  the  depletion 
region.  Yet,  in  reverse  bias,  the  voltage  drop  across  the  depletion  region  is 
large,  giving  the  impression  the  depletion  region  is  acting  as  an  insulator. 
Similar  analogies  can  be  demonstrated  when  carrier  transport  is  explained  in 
terms  of  current  flow  through  a  combination  of  insulators  and  semiconductors. 
These  analogies  will  become  clearer  upon  consideration  of  single-  and  double¬ 
carrier  currents  in  insulators  and  semiconductors. 

SPACE-CHARGE-LIMITED  SINGLE-CARRIER  CURRENTS  IN  INSULATORS 

Four  situations  are  examined  concerning  single-carrier  currents  in 
insulators: 

1.  Perfect  trap-free  insulator 

2.  Trap-free  insulator  with  thermal  free  carriers 

3.  -Trap-filled  limit  problem 

4.  Full  single-level  theory 

These  treatments  are  attributable  to  M.A.  Lampert  and  P.  Mark  (Ref.  2)..  The 
reader  can  follow  the  derivations  of  these  treatments  in  greater  depth  in 
Reference  2. 

The  basic  equations  used  for  these  various  situations  are: 

J  =  enn^E  efippE  ( 1 ) 

which  determines  the  drift  current,  but  neglects  the  diffusions  current,  and 
Poisson's  equation 

PERFECT  TRAP- FREE  INSULATOR- 

No  thermal  free  carriers  or  trap  levels  are  allowed  for  this  situation. 

The  insulator  in  question  possesses  ohmic  contacts  to  electrons  only.  Poisson's 
equation  reduces  to- 
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and  j  =en„nE  combining  with  =  o)  =  o  we  get 


/  2J  L 
E{z\  =  1  —  I  xi 


(4) 


for  I  =  i  we  get 


J  = 


(5) 


known  as  the  trap-free  square-law  relation. 


For  the  above  situation,  only  electrons  carry  current,  resulting  in 
single-carrier  current  conduction. 


TRAP- FREE  INSULATOR  WITH  THERMAL  FREE  CARRIERS 


Here  the  problem  is  complicated  by  allowing  a  finite  number  of  thermal 
carriers  to  exist  without  trapping.  Poisson's  equation  becomes 


ke  dE 


—  —  =  n  -  no 

e  ax 


(6) 


where  n„  is  some  finite  number  of  thermally  generated  electrons.  The  thermally 
generated  holes  do  not  contribute  to  current  conduction  or  space  charge  build¬ 
up.  Substituting  for  E,  and  solving  for  J  we  get 


keJ  d  (  \  \  . 


no 


A  Tik  Iv 
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kiJ  (  -1  dn{x)\  _  .  . 


^  =  (n{z)  -  no)^dx 


<in{x)  _  no  \ 

(n(*))"  V  n(x))keJ  * 


» »n  ^  j  I 


multiply  both  sides  by  nl 


n(j|  _  g'n,^/. 


let  «  = anc  t  =  and  substituting  wl  =  - - 

n(x)  /ctJ  .v.y-  ^ 

in  the  last  equation  with  w=  f udw 
Integrating,  the  solution  is  found  to  be 

with  n  =  0  and  u;  =  0  for  £(0)  =  0  w  =  —u  —  /n(l  -  u) 

Thus, "for  all  finite  w,  u  is  always  less  than  unity  so  there  will  always  be  a 
finite  amount  of  injected  space  charge  in  a  solid  of  finite  dimensions.  Sub¬ 
stituting  for  u  and  w  in  the  equation  gives 


j  =  J_ 

ke  wi 


tnoL^  vi 


Numerically  tabulating  w  versus  u  and  v  versus  u,  gives  Table  1, 
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TABLE  1.  Numerical  tabulation  of  w,  u,  and  v. 


3.61  0.99  0.277  3.125  0.99  0.239 
1.40  0.90  0.714  0.997  0.90  0.508 
0.809  0.80  1.23  0.489  0.80  0.747 
0.503  0.70  1.988  0.258  0.70  1.01 
0.193  0.50  5.18  0.068  0.50  1.82 
0.056  0.30  17.85  0.011  0.30  3.50 
0.023  0.20  43.47  0.003  0.20  5.67 
0.005  0.10  200.0  0.0003  0.10  12.0 


The  total  current  is  such  that  J  is  proportional  to  1/w  as  u  varies.  The 

9  2 

value  V  is  proportional  to  v/w'^  as  u  varies.  A  graph  of  1/w  versus  v/w  (Fig.  1) 

shows  an  ohmic  region  and  a  square  law  region  before  space  charge  limited  current 

f 1 ows . 

A  final  quantity  of  interest  here  is  the  distribution  of  charge  in  the 
material.  The  quantity 

n  ( *- ) 

•TTm 

represents  the  ratio  of  the  injected  charge  to  the  charge  at  the  end  of  the 
device.  Equation  8  is  graphed  in  Reference  2  and  is  shown  here -(Fig.  2). 

Equation  8  is  exact  at  the  extremes  of  w. ,  but  is  only  approximate  in  between 

a 

those  extremes  if  the  variables  are  expanded  in  terms  of  series. 

TRAP-FILLED  LIMIT  CASE 

In  this  case,  the  injection  level  is  such  that  the  available  traps  are  all 
filled.  Poisson's  equation  takes  the  form 


Here,  n^  refers  to  a  trapped  charge,  such  as  that  found  in  a  depletion  region 
of  a  device.  Dimensionless  variables  similar  to  the  last  case  are  used  to 
solve  the  following  analytic  equations: 


.vj'* 


i/l  »f 
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n(x)-no 
n(  L)-no 


J!!a=5.0 


10^  TN 


,V=l-05 


w  =  0.20 


0  0.2  0.4  0.6  0.8  1.0 

J. 

L 


Figure  2.  A  plot  of  charge  versus  length.  Note  the  build-up  of  spare  charge 
at  Increasing  levels  of  carrier  Injection. 
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U« 


n,  _  enttirE[x) 

i(i)  J 


wt 


e^n^fA^x 

ekJ 


■  ekJ^ 


Where  n,  is  the  electrons  mobility  (independent  of  e).  The  differential 
equation  in  E  becomes 

Ufdui 

— —  =  dw, 

1  +  Uf 

the  solution  is  wt  =  ut  -  ln(i  t  ut) 

Where  ut  =  o  at  w,  =  o  corresponds  to  E(0)  =  0. 

The  potential  is 

.  vt-  J  utdwt  =  -Uf  -  tit  +  ln(l  +  Ui) 

Jo  2 


resulting  in 


and 


J  _  e^ntftL  1 
tk  Wta 


entL^  Via 
k€  u;?a 


(11) 


where  the  subscript  a  denotes  x  =  L.  Solving  for  w  and  v  numerically  through 
u,  gives  the  values  shown  in  Table  2. 


TABLE  2. 

Numerical 

solutions  to  w 

and  V  through 

known  u. 

u 

w 

1/w 

V 

v/w 

inn 

4.9x10 

20133 

3.3x10 

134.11 

0.05  . 

1.2x10 

826.5 

4.0x10 

27.44 

0.1 

4.6x10 

213.2 

3.1x10 

14.1 

0.5 

9.4x10 

10.58 

3.04x10 

3.40 

0.9 

0.258 

0.387 

1.47x10 

2.2 

0.999 

0.306 

3.26 

1.92x10 

2.05 

0.99999 

0.306 

3.25 

1.93x10 

2.05 

1 

0.306 

3.25 

0.193 

2.05 

5 

3.2' 

0.311 

9.29 

0.902 

10 

7.6 

0.131 

42.39 

0.733 

50 

46 

0.021 

1203 

0.567 

100 

95 

0.01 

4904 

0.539 

9 
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The  hallmark  of  the  trap-filled  limit  law  is  a  large  change  in  current 
over  a  small  change  in  voltage  (Fig.  3). 

FULL  RECOMBINATION  LEVEL  THEORY  USING  THE  REGIONAL  APPROXIMATION  METHOD 

An  aspect  of  space-charge-limited  current  flow  that  can  be  used  to  explain 
the  physical  conditions  in  a  semiconductor  comes  from  the  monotonically 
decreasing  space  charge  that  naturally  occurs  in  an  insulator  (Fig.  4). 

If  the  injection  level  is  not  high,  there  will  be  a  plane  of  Xj  where 
(x)  =  no.  This  plane  will  be  a  function  of  the  level  of  injected  current 
and  will  define  two  regions  in  the  length  of  the  material.  One  region  between 
0  and  Xj,  where  ni  >  no  and  a  second  region  between  x^  and  L  where  n„  >  n, . 
Therefore,  the  simplifying  approximation  can  be  made  that  in  Region  I  (defined 
as  0  <  I  <  xi  )  n(x)  »  nj(a:)  and  in  Region  II  (defined  as  xi<x<L)  n(r) »  n,, . 
Thus,  this  method  is  called  the  Regional  Approximation.  Poisson's  equation  and 
appropriate  boundary  conditions  then  allow,  the  formulation  of  solutions  in  each 
region  showing  when  transitions  to  various  conduction  regimes  occur. 

This  method  can  be  applied  to  the  situation  of  a  set  of  traps  lying  in  an 
insulator  well  above  the  fermi  level  by  several  kT.  This  situation  can 
be  analyzed  by  constructing  a  schematic  energy  band  diagram  (Fig.  5)  that  is 
analogous  to  the  diagram  shown  in  Fig.  4. 

This  analysis  closely  follows  that  of  Lamport  and  Mark  (Ref.  2)  who  have 
used  these  methods  extensively  in  analyzing  space-cliarge-1 imited  and  double- 
injection  currents. 

Region  I  (defined  as  o<i<zi  )  is  expressed  by  the  following  expressions 
J  —  e/tnE 


n(xi)  =  Nt 
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Figure  3.  Dimensionle 
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concentration  decreases  across  the  device  F(x)  approaches 
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Region  II  (defined  as  xi<x<x2) 


where 


n(i)  =  Ncexp 


F{x)  -  E,(x) 


l  +  l/,(iNrn{i)) 

,V  -  *V,«P- 

Here,  g  stands  for  the  statistical  weight  for  one  electronic  level  to  hold  one 
electron.'  Under  certain  circumstances,  g  then  takes  on  such  attributes  as  the 
thermal  ionization  rate  per  unit  volume.  Under  present  circumstances  the 
variable  will  not  be  identified  to  a  particular  mechanism  so  as  to  leave  the 
equations  general . 

Region  III  (defined  as  <  i <  23 )  gives 


~k€  /  dE  N  n 
t  \dx  )  6 


n{x3)  =  "0 

Region  IV  (defined  as  x^kxkL)  gives 
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It  should  be  noted  that  Lampert  and  Mark  fRef.  2)  assume  that  «  «  ffnt/gNt 
until  nt  =  iv« .  That  is,  for  n<y/g,  n  nt/gjv,  and  for  n>N/g,  nt  =  N,. 

By  using  the  dimensionless  variables  used  in  the  trap-free  insulator  with 
thermal  carriers,  repeated  here  for  clarity,  the  above  Possion's  equations  can 
be  solved. 


_  no  _  enonE(x) 
“  n(i)  ~  J 


w  =  - 


V  - 


kiJ 


k€J^ 
For  Region  I : 


^  =5  i  — .  tidu  =  dw 
dw  u 


2“ 


=  r  udw=  r 

Jo  Jo 


u^dxi  =  -tt* 


«i  =  tt(wi)  =  1 


For  Region  II: 


du  Nt  no  ,  .  no 

—  - -  —  dw  — u  =  w 

dw  no  Nf 


U  =  Ui  =  1 


1 

-1  =  2 


For  Region  III: 


du  1  a  J  J 

—  =  - - -  9udu  =dw-*8  —  =  w  —  6  —  =  y 

dw  9u  23 


For  Region  IV: 


du 

—  =0-.u=  1 
dw 
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Associated  with  these  regions  are  critical  current  densities  defined  so 
that  their  appropriate  borders  (xj,  x->,  and  X3)  equal  L.  Thus,  the  critical 
currents  are  defined  as  X3[j.,ri)  =  L  x2{Jct2)  =  L  ,  xsiJcrs)  =  L .  J  <  (Ohm's  Law 
regime) . 


Examination  of  the  boundaries  show  that  for  X3(j}  <  £  all  four  regions  are 

X  V  ^ 

in  the  insulator.  Yet,  calculations  show  that  —  = <  1  and 


III  VIIW  lll,^VIIMWWI  •  I  I  Irt  f  WIIW  WIIWVV  VilUW  —  -  ..  ^  -  uiisa 

X2  ,  ^2  2 

—  =  2(pno/.V)- <  1  ,  since  SO  Regions  I  and  II  contribute  negligibly  in 


this  regime.  Thus,  Regions  III  and  IV  can  be  expanded  to  the  cathode.  The 
total  potential  drop  is  now 


V  =  Via  +  ^  =  Vui  +  (to  -  toi) 

S  ’  (i)  -  5  (i)’ 


26  -  36 


+  tOa  =  to„  -  - 

o 


and  is  tr.e  current  voltage  characteristic.  This  is  Ohm's  Law  with  a  small 
correction  term.  * 


This  area  terminates  at  a  critical  current  and  voltage  determined  by 


X3{Jerl)  =  L  : 


j  _  2e^nlnL 
6tK 


AenoL^ 


from  the  before  defined  dimensionless  variables  w  and  v. 


The  presence  of  trapping  delays  the  departure  from  Ohm's  Law  till  a 
critical  current  and  voltage  are  reached  that  is  1/6  times  the  trap  free  criti¬ 
cal  current  and  voltage. 


J„i  <  J  <  Jct7  (Shallow-Trap  Square-Law  Region) 
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Because  xsiJ)  >  l  >  <3nly  Regions  I,  II,  and  III  are  considered.  Since  the 

approximation  —«-■«:  i  is  still  valid.  Region  I  is  again  ignored.  Region  II 
^2  2 

now  extends  to  the  cathode  and  Region  III  to  the  anode.  Lampert  and  Mark 
(Ref.  2)  arrive  at 


1.,  n?,g-  1  y  _  8  /  \ 

r.  V  2NtNw.j 


(14) 


The  shallow-trap  square-low  current-voltage  characteristic  terminates  at 

X2(Jrr2\  =  I  .  SO  that 

t-nlnLN,N 
ktnlg 

eL-Nt 


Jlt2  — 


V.r2  = 


2lct 


(15) 


where  J^r2<  J  <  Jcrz  (Trap-Filled  Limit  Region). 


Since  Z2(J)  >  L ,  only  Regions  I  and  II  are  considered.  This  is  identical- 
to  the  trap-filled  limit  problem  already  observed;  i.e., 

-  -  ^  _  JVt_  _  no  _  no  /  no'  \ 

2no  u)„JVt  WaNt  \ujaNt/ 


Ji-r3  — 


2e'^N^nL 


V'.ra  = 


kt 


Zkt 


where- <  y  (Trap-Free  Square  Law  Region) 


Now  only  Region  I  exists,  which  corresponds  to  the  past  case  on  the  trap- 
free  square-law  formula. 


Figure  6  shows  the  graph  of  current  versus  voltage  and  the  associated 
regions  at  current  conduction.  The  regional  approximation  method  thus  provides 
a  powerful  tool  to  examine  charge  injected  into  a  material. 
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TWO-CARRIER  CURRENTS  IN  INSULATING  MATERIALS 

In  the  previous  discussion,  single-carrier  injected  currents  were  limited 
by  the  formation  of  a  space-charge  region.  If  the  injecting  contacts  are 
allowed  to  be  ohmic  for  both  electrons  and  holes,  space-charge  formation  is 
eliminated.  Both  holes  and  electrons  contribute  equally  to  current  flow  and 
the  phenomenon  called  double  injection  occurs.  Since  space-charge  formation 
can  not  occur  to  limit  current  flow,  a  new  limiting  mechanism  appears  to  pro¬ 
vide  a  limiting  regime  of  current  flow;  i.e.,  recombination  through  localized 
recombination  kinetics. 

Recombination  can  occur  via  recombination  centers  that  alternate  between 
capturing  electrons  and  holes.  Specialized  recombination  centers  (electron 
and  hole  traps)  may  capture  only  one  type  of  carrier  to  become  neutralized. 
Such  sices  will  be  limited  recombination  centers  whose  availability  is  depen¬ 
dent  upon  their  relative  position  in  the  forbidden  gap.  Direct  recombination 
between  valence  and  conduction  band  sites  will  also  provide' a  funnel  for 
carrier  recombination.  This  last  mechanism  will  be  directly  dependent  on  the 
dynamics  of  carrier  injection.  Finally,  carriers  that  do  not  recombine 
through  these  various  mechanisms  can,  under  low-level  injection  conditions, 
recombine  directly.  The  following  discussion  of  recombination  takes  the 
simplest  approach,  that  of  recombination  through  recombination  sites  deliber¬ 
ately  placed  in  the  material.  The  injection  level  of  holes  and  electrons  will 
also. be  assumed  to  be  high  enough  to  outnumber  the  available  recombination 
sites,  and  approximately  neutralize  each  other  at  the  anode  and  cathode  upon 
injection  into  the  material  being  examined. 

This  situation  is  referred  to  in  the  literature  as  plasma  injection 
(Ref.  2,  p.  210).  Discussions  in  AFWL-TR-85-116  (a  companion  report)  deal 
more  explicitly  with  the  plasma  aspects  of  situation.  For  now,  the  expression 
(plasma  injection)  will  be  used  as  described  here. 

CONSTANT-LIFETIME  PLASMA 

Assume  that  a  bar  of  insulating  material,  of  length  L,  has  a  hole 
injecting  anode  at  x  =  0,  and  an  electron  injecting  cathode  at  x  =  L.  The 
insulator  is  characterized  as  possessing  a  set  of  recombination  centers 
located  several  kT  above  the'fermi  level  at  room  temperature.  The  hole 
capture  cross  section  for  this  recombination  level  is  much  larger  than  the 
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electron  capture  cross  section.  Thus  the  electron  lifetime  is  constant  for  any 
injection  level.  A  similar  situation  is  expressible  for  holes  with  a  recom¬ 
bination  level  below  the  fermi  level.  The  present  discussion  will  only  look 
at  the  conditions  for  constant  electron  lifetime.  The  current  flow  equations 


Jp  =  epiApE 


Jn  =  enunE 


Poisson's  equation  is 


/  te  '\  dE 

(  —  —  =  p-  n 

\  t  j  dx 


The  conservation  equations  for  the  steady  state  are, 


-(OS-(s)-'— 


Here,-  t„  and  tp  are  the  electron  and  hole  lifetimes,  respectively.  The  assump¬ 
tion  here  neglects  changes  in  occupancy  of  recombination  centers  due  to 
recombination,  and  hence,  constant  lifetime  for  electrons  and  holes  (i.e.,  the 
injection  level  is  quite  high,  far  outnumbering  the  recombination  levels  in 
the  material).  The  boundary  conditions  are  the  same  as  before. 


Injected— Plasma  Regime 


in  this  regime  nmp  everywhere  and  The  subscripts  on  the  life¬ 
times  are  dropped  so  that  the  common  lifetime  is  r.  A  further  assumption  is 
made  that  r  is  independent  of  injection  lavel  or  position  in  the  material. 


The  total  current  density  J  is  a  constant. 


J  =  Jn  +  Jp  =  e{b+  l]^^nE  ;  b  = 
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Multiplying  Equation  4  by  b  and  adding  to  Equation  3  gives,  using  Equations 
2  and  5, 


E  (  —  ]  Ie(£]e]  = 

\dz )  ^  \dx /  ^  ktiipHnf 


(21) 


Elimination  of  (p-n)  is  attained  by  using  the  Poisson  equation.  Applying 
the  following  dimensionless  variables. 


w  =  _A£?_ 
Jbi^pr- 


_  _ 


ke 


ktE 


V  = 


dy  e(b+l)fipTn  Jr 


JH^fMpT^ 


(22) 


where  v  is  the  applied  potential  V(z)  =  f  £4*.  Substitution  of  Equation  22  into 

Jo 

Equation  21  yields- 


(£)<“(£)“> =0“-' 


\dw/  dy 


(23) 


with  the  solution 


(24) 


where  y  =  0atw  =  Wa  =  0  corresponds  to  z  =  O;  y  =  y,  =  y(uv)  is  the  value  of  y  at 
X  *  L.  Subscripts  a  and  c  refer  to  the  anode  (x  =  0)  and  cathode  (x  =  L). 


Integrating  Equation  24  with  respect  to  y  yields, 

^  “dy  =  j(yey*  -  ^y^) 


A  « 

.Wc  =  — t/z 
12^' 


(25) 


Thus,  the  potential  yields. 
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where  =  — «-  . 

120^ 

In  terms  of  the  dimensionless  variables. 


(26) 


(27) 


Equation  27  shows  that  J/L  is  a  function  of  Since  ^  =  6/5yc,  it  is  seen 
by  inspection  that 

We  ■  18  U?y 


or  that 


J 


—  (k€TftpHn) 


(28) 


where  r,  and  r„  are  the  average  transit  times  for  holes  and  electrons,  respec¬ 
tively.  As  a  reminder,  it  should  be  noted  that  this  solution  for  the  injected 
plasma  is  consistent  as  long  as  the  transit  times  for  carriers  are  much  smaller 
than  the  lifetime.  Transit  times  that  approach  or  are  greater  than  the  life¬ 
time  result  in  recombination  of  carriers  before  they  reach  the  appropriate 
contact.  This  case  then  results  in  space-charge  limited  currents.  This  situ¬ 
ation  can  be  directly  related  to  the  voltage  applied  to  the  material,  since  the 
above  situation  occurs  at  relatively  high  voltages.  It  appears  that  a  different 
regime  occurs  at  low  voltages. 

TRANSITION  FROM  A  ONE-CARRIER  SCL  CURRENT  TO  THE  INJECTED  PLASMA  REGIME 

Under  the  initial  conditions  for  the  plasma  injection  regime,  it  was 
assumed  that  the  injection  level  produced  more  than  enough  carriers  to  compen¬ 
sate  for  the  available  recombination  sites.  With  n  p  under  these  conditions, 
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the  equations  showed  that  a  current  density  versus  the  voltage  cubed  occurred. 
If  the  initial  conditions  are  relaxed  so  that  n  7^  p  anywhere  in  the  insulator, 
a  different  situation  arises.-  Again, 

J  =  Jn  +  Jp  =  eunnE  +  eUppE  =  eonatant  ^  29 ) 

The  values  n  and  p  are  separately  determined  in  terms  of  E  and  dE/dx.  There¬ 
fore, 


(6+  l)n=  — 
«/ipf 


kt  d 
e  dx 


(h+  l)p  = 


J 

</ip<  e  dx 


with  r  =  n/r„,  multiply  by  b,  and  adding  to  Equation  18  yields. 


~d(  dE  y ,  de  J 
ki(i„TnE- {E— }k€E-  +  —  =  0 


Using  the  following  dimensionless  variables: 


W  5= 


k€X 


dw  (E 
dy  bJT„ 


Substitution  of  Equation  33  into  Equation  32  arrives  at 


d  ,  d  y  d 


cPu  du 

dw^  dy 


therefore. 


U 


=  y  -  Vc  { 


(exp  y  -  1) 
(expy,  -  1) 


Where  y  =  0  at  w  =  w^  =  0  for  x  =  0;y  =  y^.  *  y(wc)  for  x  =  L. 


(31) 


(32) 


(33) 


(34) 


(35) 
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Integration  of  Equation  35  gives 


(36) 


w, 


-  «.(V,)  -  JM  -  (i  - 


(37) 


The  potential  is: 

f - - — -)  (y^-rS(expy-  1)  -2(expy  -  1)  -  2yexpy)  ^y^ 

Vexpy,  -  1/  '  3 


(38) 


(39) 


Vc  =  v(y<r)  =  jVc(^  "  3\/e)y*  +  (2  -  Mc)(l  -  A/c)yc 


where 


jv/  = — t — -  This  yields, 
exp  y^  -  1 


V  =  V(L)  = 


(40) 


(41) 


fy3  -  ( i  -  M.,]y-  +(2-  Mc)(l  -  M,.)y. 
iVr  -  (1  -  A/c)y?  (1  -  ‘^r-)-y? 


(42) 


For  Vr  >  i  (opposite  of  the  plasma  regime  case  where  yr  <  i  )  u;;:«y?/4 

and  V,.  so  that 

3 


Ve  _  * 
w?  3yc 


(43) 


Therefore, 

,  9fc<M„V2  ‘  (44) 

8612 


A  one  carrier  type  square  law  relationship  holds  as  long  as  remains  rela¬ 
tively  large.  As  the  injection  level  increases  and  y^  becomes  very  small,  a 
cubed  law  takes  hold  (Fig,  7). 
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Figure  7.  Transition  from  the  one  carrier  square  law  to  the  double  injection 
cube  law  in  an  insulator. 

Physically,  the  mathematics  supports  the  idea  that  (under  the  present 
assumptions)  low  field  conditions  promote  the  formation  of  a  space-charge- 
limited  current  due  to  the  short  lifetime  for  carriers.  As  the  injection  level 
increases,  enough  carriers  exist  in  the  material  to  overcome  the  ever  present 
recombination  levels  to  allow  double  injection. 

If  the  recombination  centers  are  removed  from  an  insulator,  injected 
electrons  and  holes  will  have  to  recombine  directly.  This  type  of  plasma  is 
is  known  as  a  bimolecular-recombination  olasma.  The  equations  for  these 
plasma  are  as  before, 


Jp  =  epUpE 

Jn  = 


except  r  is  now  defined  as  <  wtr  >  np  .  Where  ^  is  the  magnitude  of  the  rela¬ 
tive  velocity  between  an  electron  and  hole,  is  the  recombi nation  .cross 
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section  at  some  velocity,  and  <>  denotes  the  average  of  the 
Since  the  Injected  plasma  are  characterized  by  p  w  n,  it  Is 
p  by  n  everywhere  In  the  equations,  except  where  p-n  occurs 
following  equations  govern  Recombination: 

and 

J  =  Jn  +  Jp  e{b  +  l)fipnE  =  constant 

.iift 

b  =  — 


Using  Poisson's  equation,  the  current  equation  yields 


r-  ^  ^  1 

ix2  ■  U2«2i(6  +  l)Ai3  ^ 


Where 


f*R  = 


k(  <  WTji  > 

~2s 


Using  the  following  dimensionless  variables: 

w  =  ^ 

yi«Jk«{6(6  +  l)M3}i* 

2enR  _ 


u  = 


(6  +  l)(ipL^n  JL^ 


Substituting  Equation  33  into  Equation  32  yields 


[  d^u 

<fu»2 


enclosed  quantity, 
possible  to  replace 
Initially,  the 

(45) 

(46) 


(47) 


(48) 

(49) 


Lampert  and  Mark  (Ref.  2)  then  employ  the  dummy  variable  S  so  that 
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(50) 


So  S  *  -•  at  W  =  0,  i  =  +®  at  w  =  Wc  and  u  vanishes  at  these  limits, 
dimensionless  variable  v  is  now  defined  as. 


The 


V  = 


(2)r)i(*:e)i{i(6+ 


ZJiLhi 


then 


(x)  =  J  Edx  at  J  udw  ot  J  u  da  <x  j  u^ds 


V=V(5)=  I  dSaxp(^-p^^ 

—  oo 


finally, 


which  is  a  case  of  the  double- inject ion  square- law  regime. 
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(52) 


(53) 


(54) 


TWO-CARRIER  CURRENTS  IN  SEMICONDUCTORS 


Previously,  it  was  noted  that  excessive  numbers  of  thermal  free  carriers 
was  of  no  concern  in  insulators.  This  assumption  is  completely  withdrawn  when 
considering  semiconductors.  All  basic  semiconductor  theory  is  based  on  the 
concept  of  thermally  excited  carriers  playing  an  important  part  in  all  current 
flow  problems.  Other  important  variables  that  affect  carrier  dynamics  are 
temperature  and  carrier  injection  level. 

The  presence  of  thermally  produced  carriers  prevent  the  build-up  of  any 
substantial  space  charge  so,  for  the  analysis  performed  here,  local  neutrality 
is  assumed  to  hold. 

Under  the  circumstances  for  a  n-type  semiconductor  bar  with  a  hole 
injecting  contact  at  x  =  0  and  an  electron  injecting  contact  at  x  =  L,  the 
current  flow  equations  become 


e 


Dp^VriAp  (55) 

Jn  =  ^ ^  j 


(56) 


Poisson's  equation  becomes 


9 


(p-  po)  -  (n  -  no) 


(57) 


where  the  subscripted  carrier  concentrations  refer  to  thermal  free  carriers 
only.  Particle  conservatidn  equations  are 


e  dx 


d  .  d^n 


P-Po 

T 


(58) 
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and 


c  dx 


P-Po 

T 


(59) 


If  local  neutrality  does  not  hold,  then  r  corresponds  to  the  carrier  whose 
lifetime  is  most  constant  with  injection  level.  Also 

J  =  +  Jp  =  ejip(p+  6n)£-f  VyMp  ”  p)  =  constant  (60) 

Multiplying  Equation  60  by  b  and  adding  Equation  58  gives 

(61) 

which  can- be  rewritten  as 


eVr 


(£E\  _  _c_£  (  ^\ 

\ dx^  )  ktdx\  dx) 


+  (no-po)-^  +  2V^r7:2 


d^n 

dx^ 


(6  +  l)(n  -  no) 
t^nT 


(62) 


The  first  term  on  the  left-hand  side  (LHS)  of  Equation  62  can  be  dropped 
because  of  its  small  size  compared  to  the  diffusion  term  (the  last  term  on 
the  LHS  of  Equation  62). 

Under  high  level  injection  conditions  n  =  p  >>  Oq.  For  a  constant  life¬ 
time  {neglecting  diffusion  currents),  the  current  density  is 


J  =  tHp(b  +  l)n£  =  constant 


(63) 


and 


e  d 
kt  dx 


.dE 


<Pn 


+  (no-p«)^  +  2V-T^ 


(b  +  l)n 
t*nT 


(64) 


The  first  term  on  the  LHS  of  Equation  64  can  be  dropped  due  to  its  small  size 
in  a  semiconductor.  The  boundary  conditions  are  J„(o)  =  Oand  d„(i)  =  o  and  give 


and 


'll?) 


-eppVr 
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n(I)  = 


JLg 

2efinVT 


Where 


(65) 


( 


2V'rMnr\ 

6  +  1  J 


i 


the  ambipolar  diffusion  constant. 

Application  of  the  Regional  Approximation  method  to  this  problem  begins 
with  Fig.  8.  The  second  term  of  Equation  64  approximates  the  solution  in 
Region  I  and  the  third  term  of  Equation  64  approximates  the  solution  in  Region 
II.  Call  the  second  term  "Term  1,"  for  Region  I,  and  the  third  term  "Term  Z," 
for  Region  II.  That  is,  in  Region  I,  Term  1  dominates  the  equation.  The 
solution  here  is  a  decaying  exponential  for  n(x)  and  a  growing  exponential  for 
I*  E(x).  This  implies  that  Term  2  in  Region  II  is  doing  the  opposite  of  Term  1 

K  in  Region  I. 


In  Region  II,  Term  2  dominates,  giving  the  solution  of  a  monotonical ly 
increasing  E  with  increasing  x.  Since  the  solution  to  Regions  I  and  III  are 
the. same,  there  is  a  plane  where  E  is  maximized  in  Region  II  and  decreases 
into  Region  III  to  satisfy  the  boundary  conditions.  Lamport  and  Mark  (Ref.  2) 
have  solved  this  problem  and  found  that  for  increasing  voltage,  a  critical 
voltage  is  reached  where  Region  II  disappears  and  Regions  I  and  III  join. 

That  voltage  is 


v„  «  Vt 


xbi  L 
6  +  l**P2La 


(66) 


x  =  o 


X  =  L 


Figure  8.  Illustration  of  the  insulator  under  consideration  with  depiction  of 

the  general  regions  detailed  as  regions  of  differing  charge  injection. 
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Beyond  this  voltage  the  current  increases  without  further  increase  in  voltage. 
This  conclusion  shows  that  the  semiconductor  square  law  holds,  but  that  the 
length  is  replaced  by  an  effective  spacing  Leff  so  that. 


•/  j  j  e(no  -  Po)MnMpT 


te) 


where  -  xi 


For  an  insulator  injected  plasma  at  high  injection  levels  n  =  p  and 
carrier  lifetime  is  constant.  The  current  equations  are 


J  =  l)nE  =  constant 


_j_d_  (^\  ^  ^ 

k(  dx  \  dx  ’  ^  dx-  fi„r 


As  in  Fig.  8,  Regions  I  and  III  are  dominated  by  Term  2  while  Region  II  is 
dominated  by  Term  1.  Lampert  and  Mark  solved  this  set  of  equations  to  find  a 
similar  critical  voltage  as  sound  in  the  semiconductor  injected  plasma  problem. 


The  current  was  found  to  follow  a  cubed  law  as  found  with  insulators. 


(S) 


L,ff  =  X2  —  Xi 


For  semiconductor  devices  or  materials  with  heavy  doping  experiment  has 
found  excellent  agreement  as  shown  in  Fig.  9  (from  Ref.  6). 


Materials  with  very  light  doping,  such  as  the  intrinsic  region  of  a  PIN 
device  exhibit  a  cubed  law  relation  as- shown  in  Fig.  10  (from  Ref.  7). 


The  previous  discussion  lays  the  basis  for  studying  carrier  injection 
which  includes  effects  due  to  changes  in  occupancy  of  recombination  centers  or 
changes  in  carrier  lifetime  with  an  injection 'level . 


This  discussion  also  point  out  that  the  depletion  region  of  a  device  can 
be  modeled  as  semiconductor  or  insulator-like  regions  dependent  on  the  injec¬ 
tion  level  applied  to  the  device. 


DOUBLE  INJECTION  WITH  TRAPPING— THE  NEGATIVE  RESISTANCE  CASE 

Figures  11  and  12  detail  the  situations  that  will  be  analyzed  for  recombi¬ 
nation  centers  that  are  fully  occupied  in- an  insulator. 

In  this  situation,  the  capture  cross  section  for  holes  is  much  larger  than 
for  electrons  (since  full  acceptor  states  are  involved).  The  equations  govern¬ 
ing  current  flow  are 

J  =  eUnXiE  T  efippE  =  con6t>int  (71) 

the  local  neutrality  condition,  neglecting  space  charge  is 


^  =  P^Pn  (72) 

P.r  =  concentration  of  emoty  acceptor  states 

;.,(£)(»£)=. --^,.(i)(p£)  (73) 

r  -  ~  ^  ^ 

'f 

1 

—  =<  V<7„  >  PR 

•n 


1 

—  =  <  t/<7,,  >  n  /? 

(74) 


(75) 


PR  ^  nr  =  y 


Where  Nj^  is  the  total  number  of  acceptors. 
Equation  73  can  be  rearranged  to  yield 

£[(„_p)£]=  (±±^ 


Mp^p 


Mn  <> 


The  boundary  conditions  are  as  before,  E  =  0  at  x  =  0.  Using  the  Regional 
Approximation  method,  the  insulator  is  divided  into  two  regions:  Regions  I 
and  II.  In  Region  I,  extending  from  the  anode  to  x  =  x^,  p(x)>N[^  and  the 
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CONDUCTION  BAND 


FERMI  LEVEL 


ACCEPTOR  LEVEL 


VALENCE  BAND 


tR 

h 


Figure  11. 


Energy  band  diagram  showing  the 
relative  positions  of  the  acceptor 
and  Fermi  level  to  the  conduction 
and  valence  bands. 
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I 


II 


I  III  I  IV 

'  I  TRANSITION  I 

NUETRALITY  HO LDSh*.™'  '.^ELECTRON 

I  ,  I  region 

Hasp 


X  =  0 


p  -^p«n^(J«Jn)inR  =  NR 

I  i  i 

X,  X2  X3 


CATHODE 


Figure  12.  Schematic  diagram  of  the  regions  and  boundary  conditions 
for  the  study  of  double  injection  with  trapping. 


recombination  centers  are  mostly  unoccupied.  Region  II,  extending  from  the 
cathode  to  xi,  contains  fully  occupied  recombination  centers.  In  Region  I 
where  n>iV_  and  prvNk  in  Eguation  72,  Eguation  74  becomes 


=<  >  ffR 


(76) 


where  is  the  high  injection  level  lifetime,  while  Eguation  56  evolves  to 


J  =  ennE(a  +  l)n  -  eppNxE 


(77) 
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Substituting  (a  +  1)  in  Equation  77  into  Equation  76  yields 

[—L^dE  ^  —  dx 


with  the  solution 


E- 


In 


solving  for  the  potential  drop 


X 


MnT. 


(il 

JE 

(  y  - 

V  2 

W,Nr) 

enjs  «  )  \ 


At  the  edge  of  Region  I  where  x  =  X] 


"1  =  n(xi)  =  Nr 


4 


El  *  E{xi] 


J 

etin^R 


aJ 

eupNR 


=  a  J* 


Substituting  Equation  81  into  Equations  79  and  80  yields 


(78) 


(79) 


(80) 


(81) 


-I  =  ^„rpj{a  -  ln(l  +  a)} 


V(xx)  =  >i«TKJ"{^a2-«»+la(‘»+m  (82) 

Letting  and  correspond  to  Xj  =  L,  that  is,  n  =  at  the  ca inode,  then 
Equation  82  becomes 

j  m  a  tNpL 

"*  o-Iji(l-f-a) 


^  a(^a^-a  +  ln(lH-a)  (83) 

(a-ln(l  +  o))2 


after  application  to  Equation  78. 
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<  vOf  >  I*  Ntf 


ViH 


]0  w 


Examination  of  the  first  coefficient  in  a,  the  second  term  in  Equation  90, 
shows  that  the  coefficient  varies  from  about  3/2  for  a  «  1  to  2  for  a»  1.  In 
comparison,  the  coefficient  in  a  to  in  Equation  83  decreases  from  about  3/8 
at  a  «  1  to  1/9  for  a  »  1.  Therefore,  unless  a  is  very  large  >  V^,. 

Equation  90  can  then  be  simplified  further  to 


v»v„ 


(91) 


This  equation  demonstrates  that  as  J  varies  from  zero  to  Jm,  V  decreases 
from  to  Vm  where  V^h  corresponds  to  when  the  average  hole  transit  time  is 
equal  to  twice  the  hole  low  level  injection  lifetime  (Eq.  90).  Although  this 
interesting  result  still  requires  a  physical  explanation,  we  will  resolve  the 
problem  with  space  charge  for  further  insight. 

By  including  space  charge  in  this  analysis,  the  neutrality  equation  used 
in  the  previous  analysis 


n  =  p  +  pjt 

must  be  replaced  by 


(92) 


p+PR-n 


(93) 


in  certain  regions  of  the  material.  Figure  12  is  used  to  facilitate  this. 


The  space  charge  cannot  be  neglected  beyond  xz,  and  Poisson's  equation 
takes  the  form  of  Equation  93.  Hole  penetration  is  negligible  beyond  X3  to 
the  cathode  and  the  mobile  electron  charge  dominates  the  current  flow.  The 
electric  field  increases  in  strength  as  x  proceeds  from  the  anode  or  cathode 
and  peaks  in  the  transition  region. 


In  Region  III,  the  following  inequalities  hold,  p«  n,  n  ^  P(^,  p«  Hr  and 
nR  %  Nr.  The  differential  equation  for  pp  in  this  region  is  directly  obtained 
from  Equation  86,  i.e.. 
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Pr  \  aJ  J 
Integrating  Equation  94  gives 

-J: _ 

Ph(*)  P«{*a)  \  aJ  ) 

In  Region  IV,  the  following  inequalities  hold:  pc  Pj^c  N,  n^xN^. 


(94) 


(95) 


J  =  eu„  nE 


i  kt  \  dE 

"'-iTjs 


(96) 


(97) 


The  solution  found  previously  for  a  space-charge-limited,  trap-free  cur¬ 
rent  was 

i 


£:(x)=  (2L\ 


and 


(98) 


Mark  and  Lampert  (Ref.  2)  then  introduce  the  dimensionless  variables  of 
Equations  99  and  100 


I  =  w 


E  =  E-^-^u 
V  =  V*-*v 


(99) 


£+■*•  = 


2aJ 


kaWR  <  I/O-,,  > 


V  +  +  =  = 


k^t^SR  <  vop 


(100) 
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Region  I  can  be  solved  using  the  variables  of  Equation  99  and  100  with 
Equations  78  and  79. 

-  udu  <  VOn  > 


u  -t  A 

a - ai 

<  vCf,  > 

J 

ke  <  vOj,  > 

2eanp 

<  van  > 

!  <  uffp  > 

a  = 

ffjL 

l*n 

(101) 


equivalent  to  Equation  93.  The  solution  is,  like  Equation  94 


U  — 


_ In  I  \  = 


(102) 


Solving  with  the  boundary  conditions,  similar  to  those  used  to  solve 
Equation  94  results  in 


“'1  *  {a-ln(l  +  a)}' 
Ui  =  u(u;i)  =  aA 


(103) 


The  dimensionless  voltage  across  Region  I  is 

‘'(»"i)=  {^a^-a  +  ln(l  +  a)} 

In  Region  II,  Equation  86  becomes,  in  dimensionless  variables. 


(104) 


du 

dw 


(105) 


with  the  solution  u  -  Uj  =  w  -  w^,  from  Equation  103. 
Also,  the  plane  X2  in  Region  II  was  chosen  so  that 


PR. 2  *»  ’*R.2 


(106) 
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Which  corresponds  to 

(  \  1 

- 

«>3  «  "  Oi4  +  {0  -  ln(l  +  a)}’ 

(107) 

•1  ttl 

Wi.3  j  udw  =  j  udu  =  i(u|  -  «?)  =  ^  (j  -  j 

•l  «l 

(108) 

To  evaluate  Equation  95  at  x  =  x^,  fix  P^jto  be 

(109) 

since 

p.., -!»(.,) --i  (7)  (g)^ 

(110) 

£  is  used  to  fix  x,.  Equation  95,  in  terms  of  the  dimensionless  units  becomes 

or 


4(u;r  —  u»3)^  =  1  +  2(u)3  —  wt) 


(111) 


To  approximate  u(E)  in  Region  III  an  interpolation  scheme  is  used.  Let 
u  *  U2(w)  be  the  equation  for  the  tangent  to  the  curve  u  =  u(w)  in  Region  II 
at  w  “  W2.  Let  u  *  U3(w)  be  the  tangent  to  the  curve  u  =  u(w)  in  Region  IV 
at  w  »  W3.  Then 


(w)  m  U2  +  (u,  _  »  Uj  +  -  Wj) 


(112) 


and 


“3  »  “3  (w  -  *V3j 


“■''(si*  ■ 

u;  -  ti’ 

2(u;,  ~  Wa)! 
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Note  that 


du- 

dw 


=  -1 


(114) 


SO  that 


u  =  (to,  —  w)  i 


and 


du  1  . 


for  the  boundary  condition  u,  =  u(ti;<;)  =  o  at  E(L)  =  o. 
The  voltage  applied  to  Region  IV  is 


udui  *  -  u(3)i 


(115) 


By  using  Equation  100  we  have  the  relations  determining  the  current-voltage 
relationship  can  be  determined.  They  are 


2anpJ 


w,  ktNii'  <  wXp  L  ■ 


urf  Sg  <  uffp  >-  L- 


(11.6) 


Lampert  and  Mark  (Ref.  2)  solve  these  equations  in  terms  of  W  and  V  to  be 


-1  -  J.4  -  (i)  {(a  -  ln{l  a))  -  1)}  -  (v+  l)= 


(117) 


i;,  |y^+  4»-  (ii?)  (118) 

w?  [-1  -aA^  (|)  {(a  -  b(l  a))  +  1)}  +  (y  +  I)®]* 


Where  y  *  (w^.  -  W3)  *  and  V^-  =  V^^l  +  2  ^  '^2  3  '^S.c’  substitution  in 

Equation  93  then  yields 


(i:). 


k{a)b 


aA 


and 


(±)  =_1_ 
\»?/m  ®«(®) 


(119) 


di 
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where 


^  _  ke  <  v<7p  >  _  J 

2eanp  efipNaE-^-^ 


<vcp>  iiprpE++ 


h(a)  = 


a  —  In  (1  +  a) 


_  {a-In(li- g)}^ 
a{i  =  -a  +  ln(l  +  a)} 


This  allows  the  equation  (Eq.  91)  for  the  neutrality  based  characteristic 
to  be  written  as 


(120) 


The  space-charge-based  J-V  characteristics  are  compared  in  Fig.  13  to  the 
neutrality-based  J-V  characteristics  for  two  sets  of  material  parameters  in 
silicon.  The  lower  set  of  curves  correspond  to  k  =  12, 


ftp*  10*; 


<wa,>  =  9x  10“‘®cm^/s 


<  vap>  =3x  10-«cmVs 


These  values  are  for  silicon  near  liquid  nitrogen  temperature.  The  dimension¬ 
less  variables  correspond  to  a  =  1,  A  =  lO"^  and  B  =  3  X  10"**.  The  upper  set 
of  curves  correspond  to  A  =  3  X  10"**  and  B  =  10"3.  These  situations  show  that 

a  change  in  carrier  lifetimes,  r= - -  ,  directly  determines  when  the 

negative  resistance  regime  will  occur.  This  can  even  be  brought  about  by  a 
change  in  material  temperature  affecting  or  <  v<t>  . 
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Figure  13.  Current  density  versus  voltage  for 
the  case  of  double  injection  in  an 
insulator. 

Examination  of  the  four  regions  assumed  for  this  material  as  the  lower 
J-V  curve  is  followed  also  shows  some  interesting  aspects.  At  a  on  the  lower 
curve.  Region  IV  occupies  83  percent  of  the  material  and  81  percent  of  the 
applied  voltage.  Region  III  takes  the  remaining  16  percent  of  the  material  and 
18  percent  of  the  voltage.  Regions  I  and  II  take  up  less  than  one  percent  of 
the  material  and  voltage.  At  point  6,  Regions  I  and  II  take  up  32  percent  of 
the  material  and  about  two  percent  of  the  voltage.  Region  III  takes  43  percent 
of  the  material  and  72  percent  of  the  voltage.  Region  IV  takes  25  percent  of 
the  material  and  about  26  percent  of  the  voltage.  At  y.  Regions  I  and  II  pos¬ 
sess  55  percent  of  the  material  but  only  five  percent  of  the  voltage  compared 
to  Region  III:  34  percent  of  the  material  and  80  percent  of  the  voltage. 

Region  IV  11  percent  of  the  material  and  five  percent  of  the  voltage.  So  far 
Region  IV  is  continuously  shrinking  in  size.  At  6,  Region  I  and  II  occupy  99 
percent  of  the  material  and  98  percent  of  the  voltage.  Region  IV  is  negligible 
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and  Region  III  takes  the  residual.  The  material  moved  from  a  space  charge 
limited  current  flow  to  one  of  plasma  injection,  demonstrating  the  theoretical 
existence  of  a  negative  resistance  regime.  As  double  injection  occurs  an  elec¬ 
trically  unstable  situation  develops  and  the  material  moves  to  a  more  stable 
state  as  the  J-V  curve  breaks  over  into  current  controlled  negative  resistance. 

RECOrCINATION  CENTERS  PARTIALLY  OCCUPIED 

In  contrast  to  the  previous  problem,  the  recombination  centers  in  this 
case  are  allowed  to  be  only  partially  occupied.  This  can  be  accomplished  by 
placing  the  recombination  centers  near  the  Fermi  level  and  assuming  that  the 
number  of  free  carriers  are  negligible,  i.e.,  the  centers  are  primarily  defect 
centers  or  charged  traps. 

Once  again,  the  capture  cross  section  for  holes  is  assumed  to  be  much 
larger  than  the  capture  cross  section  for  electrons.  As  a  result,  many  of  the 
same  characteristics  seen  in  the  prev.ious  problem  will  occur  here. 

The  basic  difference  in  this  problem  is  the  occurrence  of  a  recombination 
barrier  to  current  flow.  In  the  previous  situation,  only  one  carrier  type  saw 
any  obstruction  to  current  flow.  Under  present  conditions  both  electrons  and 
holes  see  a  voltage  barrier  that  must  be  overcome  before  current  can  flow 
(Fig,  14)’. 


CONDUCTION  BAND 


TRAP  L*EVeT^*^  ^  •- 

- Eu 

VALENCB  BAND  ^ 


Figure  14-.  Energy  band  diagram  for  the 
case  of  injection  with 
recombination  centers  par¬ 
tially  filled. 
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Since  trapped  space  charge  exists,  the  neutrality  approximation  must  be 
dropped  and  Poisson's  equation  now  assumes  the  form 

(7)(f  <121) 

where  p«.-  is  the  empty  recombination  centers.  In  addition,  the  boundary  con 
ditions  take  the  form  E  ■  0  at  x  *  0  and  x  *  L.  The  electron  flow  dominates 
at  low  current  levels  so  that 


J  =  en„  nE  V  conatant 


(12?) 


and 


Since 


then 


n  p  n  p 

T„  Tp  <  V<T„  >  Pit  <  VCp  >  ng 


(123) 


(124) 


p  <v<rp>  ng  -  n  <van>  pR 

Multiplying  by  e  E  and  substituting  from  Equation  122  results  in 


Jr  *» 


Mn  \<VOp>  /  \ngj 


Since  PR-rng  =  yg  differentiating  ^  gives 


(125) 


(126) 


±  (f£\ 

dx\ng)  \n%)  dx  \  e  )  n%  \dx^  J 

Differentiating  Equation  126  and  substituting  Equation  127  results  in 


(127) 


djp  I <Wn>  d'^E 
di  **  I  e  <  wTp  >  n  J  )  dx^ 


(128) 


Simplifying  gives 
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This  can  be  used  in  Equation  127  to  obtain  . 

E  ( ss  *  ^  e  <vOp>  n%  „pR,o 

\dx^  )  ^  ken,, Nr  kenpNn  (130) 

This  implies  a  threshold  voltage  independent  of  current  level.  Solving  Equation 
130  via  simple  approximations  (from  the  low  current  conditions)  yields. 


e  <  v<T,,>  L*n%  ,^pR.,j 
kefipNR 


(131) 


a  voltage  barrier  that  must  be  overcome  for  the  current  to  increase  appreciably. 
Physically,  this  corresponds  to  a  situation  where  the  recombination  centers  must 
be  filled  to  a  major  extent  before  current  can  flow — a  trapped  space  charge 
must  be  overcome. 

Further  insight  can  be  acnieved  through  application  of  cne  Regional 
Approximation  method,  i.e.,  dividing  the  sample  used  into  three  sections.  The 
third  region  in  Fig.  15  corresponds  to  a  region  described  by  the  most  recent 
set  of  calculations,  where  Region  III  represents  a  potential  barrier  (i.e., 

'^th  dissipated.  Local  neutrality  is  assumed  for 

Regions  I  and  II.  In  Region  I,  n{x)>nf^^Q  and  the  electrons  are  primarily  in 
the  conduction  band;  while  in  Region  II,  n(x)  <  n|^^Q  and  the  population  of  the 
recombination  centers  is  not  changed  from  the  thermal  equilibrium  values. 

-  In  Region  I,  n(z)  >  n.>?. .  «  .Vr  .  Local  neutrality  is  detailed  by 

^  ~  P  ~  ~~  PR,ff  “  ^R.a  f  .  oo  \ 


n  =  n. 


n  =  JlL  ^ 
e  dk 


Figure  15.  Schematic  diagram  of  regions  under  consideration  with 
boundary  conditions. 
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Thus, 


feiO" 


—  =<Vffn> 


Jt  =  ei*nE(a  +  l)n  -  etXpTiR.oE 


Letting 


\  v/*p*'l iJ.O  / 


dE  =  Tdx 


the  solution  is 


r-  ,1 

-  J  ,;n  ,  .  [  =  Tx 


Using  the  boundary  conditions  E  =  0  at  x  =  0  results  in 


V(.)  =  [e,z  =  i  + 


X  =  Xi 


fii  =  n(ii)  =  nR,o  and  Ei  =  E{xi)  =  ( — - - )  =  aJo 

\^f^n^R.a  / 


Xi=  (-^)(a-ln(l  +  a)) 

Va.i  =  (^a*-a  +  ln(l  +  a)) 


(133) 


(134a) 


(134b) 


(135) 


(136) 


(137) 


(138) 


AFVIL-TR-85-115 


as  before,  solving  for  and  gives 


Jn,  =  h(a) 


(139) 


For  J  <  Jn 


(140) 


For  Region  II,  rt  r  r*  r  'R.  r  n  -  prt...  ^  n.  n/t  ss  np,,j .  The  neutrality  and  cur¬ 
rent  conditions  are 


n  =  pp  ~  pr,o  =  njf.,,  -  np 


(141) 


J-r  =  efinTiE 


(142) 


Combining  Equation  127  with  Equation  141,  and  solving  for  the  change  of  hole 
current  with  distance  yields 


aJ,,  aJ'j-  <  vffn  >  .Vp  dn 
— ^  =  - : - : — 5 - -  =  -cn  <  vo„  > 


dr  <  va„  >  ni  dx 


(143) 


Substituting  Equation  141  on  the  right-hand  side  (RHS),  and  n^  for  np  ^ 


e  <  vCp  >  nj,  ^  e  <  v<t,,  >  n^j  „ 


«("  +  PKo)  aJrNp 

JrdE  _ 


J  I 


(144) 


since  fin  =  -  ^  (from  Equation  142)  then 


-  - 

V  (ini'll .  y 

The  point  X2  is  where 


e  <  vCTp  > 


aJj-Np 


(145) 


(7)  (s),  .  . 


(146) 
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Mark  and  Lampert  (Ref.  2)  now  change  teethe  dimensionless  variables 

X  = 

E=  E^u 


V  =V^v 

(147) 

Yt  -  (iNrJt 

e  <v(Tp>  n\  ^pR_^ 

£.t  _  xNrJt 

{e  <  v(Tp>  n%  ^pR_olcefipNjt}^ 

(148) 

Vt  -  (aiVfi  Jt)^  ^  ^txt 

(149) 

For  Region  I,  Equation  134b  becomes 

**4u  _  . 

- =  Ddw 

u  +  c 

^  <V<Tp>  PR.c,'\  * 

V  ia-UpNft  ) 

(150) 

„  {  1* 
n  =<  v<r„  >  Nr  < - = - ? 

[e  <v<Tp>  n\„pR.oHp\ 

(151) 

solving  Equation  150 
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For  w  =  w(L),  Equation  153  and  154  can  be  rewritten  as 

and 

For  large  injection  levels.  u,cui,  «'c  *»  .  and 


(155) 


(156) 


Therefore , 


Jr  as  -efinfiR 


(157) 


which  is  the  semiconductor  injected  plasma  characteristic  for  n©  -  Pq  = 
Lamport  and  Mark  go  on  to  solve  Regions  Hand  III  similarly  to  obtain  the 
following  equations  for  the  current  voltage  characteristics. 


We  e  <viTp>  nji  ^PK.oL' 


ill  _  _L_I  1 

n/t.oL^  epR.o  <v(Tp>  j 


(158a) 


(158b) 


which  are  plotted  in  Fig.  16  for  the  following  parameters:  e/so  =  12,  Up  =  Up 
=  1  X  10‘*cm2/V-s,  Q  ■  Pr  0  ~  5  X  10^‘*cm"3,  <  vap>  =  lO'^cm^/s  and  <  vap  > 

=  10"^cm3/s. 

Figure  16  shows  the  general  I-V  characteristics  for  double  injection  in 
an  insulator  with  partially  filled  recombination  centers. 

It  should  be  noted  that  once  the  recombination  centers  are  filled,  the 
problem  reduces  to  that  of  the  trap-filled  limit  of  the  double-injection  regime 
Under  actual  conditions,  the  transition  from  single  carrier-injection  to  double 
carrier  injection  will  depend  on  the  rate  at  which  voltage  is  applied  to  the 
specimen.  Thus,  the  curve  may  be  more  rounded  or  abrupt  in  its  transition  to 
the  J<xV^  regime. 


Figure  16.  General  I-V  characteristics  for  double  injection  in  an  insulator. 


NOTE:  The  definitions  of  a,  s,  y  and  6  on  the  curve  are  as  follows:  a-- 
Region  I  is  1  percent  of  the  available  volume.  Region  II  is-  1  percent  of  the 
available  volume.  Region  III  is  97  percent  of  the  available  volume  and  100 
percent  of  Region  V;  s--Region  I  is  14  precent  of  the  available  volume. 

Region  II  is  7  percent  of  the  available  volume.  Region  III  is  79  percent  of 
the  available  volume  and  100  percent  of  Region  V;  y— Region  I  is  48  percent  of 
the  available  volume  and  4  percent  of  Region  V,  Region  II  is  26  percent  of  the 
available  volume  and  14  percent  of  Region  V;  and  6— Region  I  is  56  percent  and 
8  percent  of  V,  Region  II  is  29  percent  of  the  available  and  35  percent  of  V, 
Region  III  is  15  percent  of  the  available  volume  and  59  percent  of  V. 


mmimw. 
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DOUBLE  INJECTION  WITH  TRAPPING- FREE  THERMAL  CARRIERS  PRESENT 

The  same  problem  analyzed  in  the  last  section  is  now  examined  with  a 
p-type  insulator  with  thermal-free  carriers  (holes)  present.  The  same  assump¬ 
tions  apply  as  before  except  for  the  rewriting  of  Poisson's  equation. 


Jj.  =  J  =  eftnE  +  J  =  eupE  =  constant 

(159) 

(^)  (^)  ~  ” 

(160) 

(161) 

r  -  g  =  r„  =  rp  ~  g,. 

(162) 

with  n  =  0. 

II 

n 

1 

II 

(163) 

1 

—  =<U<7„>  PR 
“^n 

1 

—  =<  vop  >  nR 

1  AT 

—  =<ep>  Nv 
rR 

(164) 

Trapped  electrons  are  not  allowed  to  be  reemitted. 

Equation  159  becomes 

Jj-  =z  J  =  ennE  ^  „  j  +  tfipPoE 

(165) 

At  this  point,  the  quantity  (p  -  Po)/n  is  assumed  i 
injection  level.  This  implies  that 

constant,  independent  of  the 

n  ^ 

•  —  *  —2— 

(166) 

is  true  at  low  levels  and  high  levels  of  injection 

.  Thus, 

1 

—  =<  W(Tn  >  PR.o  =<  Vffp  >  nR,„ 

Tji 

(167) 
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and 


-2-  =  i 


(168) 


Equation  165  now  becomes 


*  -MpPo 


l-r 


n  — 1 

Mn  /p-  Po\  dE 
ltp\  n  )  dx 


(169) 


Since  free  electrons  are  trapped  quickly  without  reemission,  the  quantity 
underlined  in  Equation  169  must  be  quite  small,  resulting  in 


dx  T„  \  PR  J  r„ 


(170) 


This  equation  can  be  rewritten  in  terms  of  the  carrier  transit  time,  using 
k€ 


* 


<  V<Tn  >  PR,o 


(171) 


Thus 


\e  J  dx  r„ 

Jt  —  *A*n  —  ep,,  Pn  E  ^  constant 


(172) 

(173) 


Jhe  last  two  equations  are  similar  in  form  to  those  found  in  the  one- 
carrier  space-charge=limited  current  with  shallow  traps.  The  I-V  relations 
can  be  written  as 


J  =  epoMp 


(x) 


where 


(174) 

(175) 


or  the  transit  time  equals  the  lifetime 


“  t/  ^ 


(176) 


after  which 
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When  tp>r„,  carriers  are  trapped,  and  the  space-charge-limited  current  flow 
dominates.  For  t,,>^  double  injection  becomes  prevalent. 

Under  realistic  conditions,  the  assumptions  above  do  not  hold  as  well. 

The  lifetime  of  carriers  vary  with  injection  level  and  captured  carriers  are 
reemitted.  The  net  result  is  complicated  equations  which  must  be  solved 
numerical  ly. 

If  the  variables  are  assumed  to  change  slowly,  the  above  equations  still 
hold  generally  and  the  conclusion  that  as  the  lifetime  and  transit  time  became 
equal,  double  injection  prevails  over  single  carrier  injection  still  holds. 

This  point  is  the  one  prevailing  idea  throughout  the  mathematical  outline 
of  double-  and  single-carrier  injection.  Thus,  it  is  an  important  crossover 
point  to  understanding  the  conditions  necessary  for  current  filamentation. 
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CURRENT  FILAMENT  FORMATION 

A  notable  aspect  of  the  past  mathematical  treatment  of  various  current 
flow  problems  is  the  appearance  of  current  controlled  negative  resistance  (CCNR) 
in  I-V  characteristics. 

B.  K.  Ridley  (Ref.  8)  has  pointed  out  that  the  occurrence  of  CCNR  results 
in  the  formation  of  current  filamentation.  Therefore,  the  previous  treatment 
provides  clues  to  the  initial  conditions  under  which  CCNR  occurs.. 

To  demonstrate  the  connection  of  the  preceding  sections  to  filament  forma¬ 
tion,  Ridley's  basic  theory  is  examined  in  detail,  bringing  in  concepts  of  the 
weakly  or  strongly  ionized  solid  state  plasma  as  they  apply. 

A  bar  of  semiconductor  material  is  subjected  to  a  constant  current  and 
electric  field  and  undergoes  a  perturbation  under  constant  volume  and  pressure. 
The  second  law  of  thermodynamics  as  derived  by  Gibb  (Ref,  9)  is  as  follows 

GsV-TS  +  pV  (178) 

where 

G  =  Gibbs  free  energy 
U  =  Total  energy 
T  =  Temperature 
S  =  Entropy 
P  =  Pressure 
V  =  Volume 

Under  the  second  law,  G  is  a  minimum  for  a  system  in  equilibrium.  Thus, 
very  little  free  energy  exists  in  the  system  to  drive  instabilities  or  pertur¬ 
bations  from  equilibrium. 

The  Gibb's  free  energy  is  closely  related  to  the  helmholtz  free  energy 
allowing  the  statement  of  a  thermodynamic  identity  relating  the  two. 

dU  =  TdS  -  pdV -i- ndp  _  (179  ) 

where 

u  =  total  chemical  potential 
p  =  carrier  density 
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Substitution  of  Equation  179  into  the  differential  of  Equation  178  under 
conditions  described  above  and  constant  carrier  density  (or  current)  yields. 


d6  =  0  and 


TdS  =  dU-  ixdp 


(180) 


must  be  a  minimum  for  a  system  in  equilibrium. 


Taking  the  time  derivative  of  Equation  180  gives 


dS  du 

IT 


(ISI) 


where  the  last  term  is  for  all  species  of  carriers  (holes  and  electrons). 


The  application  of  an  electric  field  results  in  a  perturbation  of  the 
total  energy  and  the  chemical  potential  so  that 


0  =  U  -r  ^  ekPkd> 


(182) 


/Tfc  =  /*fc  +  ek<l> 


(183) 


where 


<fc  =  charge  per  particle 

<p  =  inner  potential  due  to  the  applied  field 


-The  time  derivative  for  Equation  167  is 


dU  J  •  r  , 


(184) 


where  the  total  energy  flow  is  Jy 


Ju  —  ^,ekJk<P 


(185) 


and  Jq  is  the  total  heat  flow.  Under  initial  conditions  the  heat  flow  is  equal 
to  the  electrical  energy  flow  (  ).  In  addition,  the  continuity  equation 


also  holds. 


dpk  •  ,  ,  dpk 

■jI'  *  —dtvJk  +  -T — 
■dt  •  at 


(186) 
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Substitution  into  Equation  181  yields  the  following. 


dS  dU  ^  apt 
^  dt  ~  dt  .  dt 


dpk 


r-;-  = 

dt 


f  =  +  )  -  (L  +  ^)) 

-div  ^ ekJk<^  +  gfcPt  ^1 
-  +Yi‘k4>  ^-divJk  + 

^  =  -div 


(187) 

(188) 


(189) 


(-? 


«fc*/fc  •  -  ^<fc^V  •  Jfc 


3/>fc 

dt 


(190) 


using  the  following  relationships 
E  —  -grad^ 

j  =  y  uJk 


J./  =  -  E  ^kJk<t> 


Mfc  =  -‘k<f> 

and 

^-0 

dt 

we  find 

that 

dS 

T—  =  -div 
dt 

JE  + 


(E(-' 


dpk 


k  -t-  tkdldivJk 


(191) 


(192) 
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^(Mfc  +  ek(t>)dtvJk 


=  div  ^  HkJk  -  ^  Jkgraduk  div  ^  «fc  Jfc  ^ 


=  grad{J^  -  ^  Jkt^k) 


(193) 


where  grad  =  Tgrad 


the  temperature  gradient 


therefore 


i  = 


+  Z  Jt,  grad 


■T  -  ; 

^0; 


J^u^)  grad 


(194) 


'-^hich  can  be  written  as 


dt  3t 


(195) 


Where  j,  is  the  entropy  flow  and  is  the  entropy  production  rate.  la  the 


steady  state,  is  zero  as  is  —-and  At  equilibrium  the  production  and 

dt  ■  dt  ■  it 


flow  of  entropy  is  zero  if  the  steady  state  is  quite  near  equilibrium.  Thus, 


—  should  be  as  small  as  possible  at  equilibrium  (i.e.,  minimal  amount  of  free 


energy) . 


The  concept  invoked  to  justify  the  minimization  of  entropy  production, 
namely,  the  principle  of  least  entropy  production,  is  acceptable  as  long  as  the 
definition  of  steady  state  is  not  synonynrous  with  quasi-steady  state.  For  this 
derivation,  the  time  required  to  achieve  the  steady  state  must  be  very  short 
-ompared  to  the  duration  of  the  experiment  or  observation. 


If  this  is  not  the  case  then  quasi-static  situations  prevail  resulting  in 
the  production  of  entropy  (or  release  of  free  energy)  and  ensuing  instabilities 
in  the  solid  state  plasma  formed  (i.e.,  unstable  electrical  conditions  as 
shown  in  the  companion  report). 


.Equation  194  can  be  applied,  to  a  s’emi conductor  undergoing  second  break¬ 
down  and  simplified  if  the  temperature  gradients  in  the  material  are  assumed 
to  be  extremely  small.  Equation  194  simplifies  to 
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.  ^  gfadfXieJk  ,  ^  ,  JE  ^Ukdpk 

^  -  2^  Tfc  2^  Tk  T  2^Tk  dt  (196) 

k  k  k 

Breaking  Equation  196  down  into  the  entropy  production  in  the  high  current 
region  (52)1  and  low  current  region  (s^)  upon  filament  formation,  then 

s  =ii{l-a)-i2a  +  A(i_^s*)  (I97) 

where  a  is  the  filament  cross-sectional  area  (Fig.  17). 

Two  processes  occur  upon  stabilization  of  the  current  filament:  (1)  the 
filament  will  grow  or  shrink  depending  on  the  amount  of  current  shared  in  the 
low  and  high  current  regions  and  (2)  the  entropy  production  increases  or 
decreases  (increasing  the  circuit's  entropy)  for  the  semiconductor  overall. 

If  stabilized,  the  temperature  and  chemical  potential  gradient  will  be 
zero  and  charge  conservation  will  be  observed.  From  symmetry,  the  terms  in 
Equation  196  will  cancel  except  JE/T.  Thus,  the  only  entropy  source  is  that 
due  to  joule  heating. 

T:  =  E.JAl  -  a)  -  E.2l:a 


■Jo  =  ./l(l  —  <*)  +  '^2'! 


r.-  =  EJ„  =  E{Ji{l  -  a]  ~  J2<i]  (198) 

Then,  the  steady  state  is  that  in  which  the  field  E  is  a  minimum  or 


Jo  ~  Jio 

a  = 

**70  *  JlO 

The  cross-sectional  area  of  the  filament 
cross-sectional  area  of  the  sample.  The 
upon  examination  of  Poisson's  equation 

divE  = 

dtvE  - - =  —divJ 


(199) 

is  in  the  same  units  used  for  the 
build-up  of  a  filament  is  inevitable 

(200) 

•  (201) 
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representation  of  a  current  filament. 


Integration  over  the  space-charge  layers  involved  yields. 

E.-Eo=  —{J2  -  J^) 

€«0 


(202) 


but 

then 


=  .E2  —  £.0  =  *rT(*^2  ~  *^o)  =  TT^'^2 
dJ  dJ 


•  AvdJ 


(203) 


and  a  negative  differenti-il  resistance  leads  to  the  growth  of  the  filament  to 
a  state  where  the  minimum  amount  of  entropy  is  produced  under  steady  state 
conditions. 

It  should  be  noted  that  these  derivations  assume  a  zero  phase  difference 
between  voltage  and  current,  and  a  fixed  total  current.  Dropping  these 
assumptions  adds  two  terms  and  a  phase  term  to  Equation  196.  The  added  terms 
can  result  in  driving  forces  that  can  feed  instabilities  (create  free  energy) 
and  produce  growing  oscillations  (Refs.  10,  11,  and  12)  that  influence  fila- 

4 

ment  growth  (Ref.  10). 

Couched  in  terms  of  the  previous  derivation  on  carrier  transport,  the 
achievement  of  a  transit  time  less  than  the  lifetime  creates  an  unstable 
electronic  state  that'is  mitigated  by  filament  formation.  This  situation  can 
become  stable  if  all  free  energy  (heating  or  external  drive  mechanisms)  is 
removed  from  the  device.  If  not,  instabilities  recurr  resulting  in  an  unstable 
plasma  formation  followed  by  thermal  instabilities  or  growing  oscillations. 
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APPLICATION  TO  GENERAL  CASES 
DEPLETION  REGIONS  AS  INSULATORS  AND  SEMICONDUCTORS 

As  noted  previously,  the  electrical  characteristics  of  semiconductors  are 
exploited  by  controlling  the  material's  conductivity  by  passive  (doping)  and 
active  (induced  electric  field)  means.  One  way  to  unite  passive  and  active 
modulation  of  conductivity  is  to  form  a  metallurgic  junction  between  dissimi¬ 
larly  doped  materials.  This  type  of  junction  is  generally  called  a  PN  junction 
or  diode  (if  p-type  and  n-type  materials  are  joined).  The  materials  reach 
equilibrium  after  joining  to  create  a  depletion  region.  So  called  because  of 
the  lack  of  mobile  charge  carriers. 

The  potential  developed  in  the  depletion  region  creates  an  electric  field 
that  sweeps  the  region  free  of  mobile  charges.  Injecting  a  charge  into  the 
depletion  region  via  regions  adjacent  to  it  demonstrates  a  distinctly  differ¬ 
ent  conductivity  depending  on  the  impressed  electric  field’s  polarity. 

In  forward  bias,  the  depletion  region  exhibits  a  high  conductivity  similar 
to  a  heavily  doped  semiconductor.  In  reverse  bias  the  depletion  region  exhibits 
low  conductivity  similar  to  an  insulator  (Fig.  18). 

Basically,  under  low  level  conditions,  the  conductivity. of  a  PN  junction 
can  be  modulated  over  a  wide  range  through  variation  of  the  strength  or  polarity 
of  the  applied  electric  field.  Under  high  level  conditions  the  conductivity 
is  still  modulated  but  is  dependent  on  more  than  the  electric  field.  New 
carrier  transport  mechanisms  occur;  nonlinearities  due  to  carrier  mobility 
changing  with  carrier  energy  occur,  and  variables  become  dependent  on  tempera¬ 
ture  and  device  architecture. 

For  example,  the  above  PN  junction  can  be  driven  into  reverse  bias  by  a 
strong  electric  field  and  initiate  avalanche  breakdown.  The  result  is  an 
increasing  current  capacity  at  a  relatively  constant  voltage.  The  depletion 
region  exhibits  a  constant  or  increasing  conductivity  in  response  to  the 
applied  voltage^  similar  to  a  metallic  material. 

By  taking  this  type  of  approach  to  semiconductor  transport  problems  the 
Regional  Approximation  Method  becomes  very  useful  in  solving  low  and  high 
level  device  characteristics. 
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Figure  18.  Ideal  low  voltage  I-V  characteristics  of 
a  PM  junction. 

PN  JUNCTION— FORWARD  BIAS 

A  demonstration  of  the  utility  of  the  Regional  Approximation  Method  is 
easily  performed  for  this  case.  An  abrupt  PN  diode  is  detailed  below  for 
analysis  ignoring  thermal  effects. 

P  =  N  =  1  X  10^'  carriers/cm3 

total  length  =  2  x  lO’^  cm 

cylindrical  geometry 

radius  »  5  x  10’^  cm 

u  •  r,  a  I  X  10”®  s 

«  1350  cm2/V-s 

Mr  »  480  cm2/V-s 

Built-in  potential  =  0.577  V 

Unbiased  depletion  layer  width  =  1.7  3  x  10"“*  cm 
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The  diode  is  forward  biased  and  passes  a  very  small  current  on  the  order 
of  the  saturation  current.  The  equations  used  to  model  behavior  of  a  PN 
junction  are  well  known  and  will  not  be  derived  here.  By  passing  the  low  and 
moderate  level  injection  regimes  we  now  apply  the  Regional  Approximation 
Method  to  the  high  level  injection  regime. 

Under  high  level  injection,  the  depletion  region  shrinks  to  a  fraction  of 
its  original  size.  This  follows  from  the  application  of  an  external  potential 
which  opposes  and  eventually  counteracts  the  junction's  built-in  potential. 

As  this  occurs,  the  external  potential  begins  to  fall  primarily  on  the  bulk 
regions  outside  the  depletion  region. 

As  Fig.  19  shows,  a  potential  gradient  is  impressed  on  these  regions 
forming  an  electric  field.  Charge  build-up  near  the  contacts  also  contributes 
to  an  increasing  electric  field  (Ref.  13). 

Under  low  and  moderate  injection  levels,  the  depletion  region  was  the 
major  influence  on  the  devices  characteristics.  The  situation  changes  dramat¬ 
ically  under  high  level  injection,  bringing  the  bulk  regions  directly  into  play. 
With  the  depletion  region's  influence  minimized,  diffusion  of  carriers  plays  a 
very  minor  rule  in  controlling  the  device's  characteristics.  The  electric 
field  in  the  bulk  regions  initiates  drift  and  space  charge  which  allows  a  break¬ 
up  of  the  device  into  two  major  regions  where  drift  is  dominant,  and  the  recom¬ 
bination  centers  are  full  (ionized),  as  shown  in  Fig.  20. 

The  problem  can  be  further  simplified  by  examining  one  side  of  the  diode, 
since  the  other  side  should  react  with  the-  same  transport  properties. 

Examination  of  the  P  of, Fig.  20  side  results  in  the  formation  of  four 
zones  in  Region  I  (Fig.  21). 

Electrons  are  injected  from  the  cathode  into  the  P  region  (Fig.  20)  but 
cannot  traverse  the  region  completely  before  recombining.  The  transition 
region  is  required  to  smoothly  join  Region  II  with  Region  IV.  This  situation 
is  exactly  like  the  case  of  double  injection  with  trapping  (recombination 
centers  are  full— including  space-charge  effects).  That  is,  a  JaV^  followed 
by  breakdown  to  another  JaV^  regime.  The  Fig.  22  (from  Ref.  14)  shows  the 
initial  JaV^  does  occur  under  appropriate  conditions.  The  N-side  can  be 
drawn  in  the  exact,  fashion  as  the  tn  Fig.  21  as  a  mirror  image  (except  for 
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Xg  =  METALLUGiC  JUNCTION 

Figure  19.  As  the  external  potential  increases  and  overcomes  the  built-in 

potential,  the  bulk  regions  become  involved  developing  an  electric 
drift  field.. 
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REMAINING 
DEPLETION  REGION 


0  Xq  L 


Figure  20.  Depiction  of  a  PN  junction  under  heavy  forward  bias  (not  to  scale). 


Figure  21.  Separation  of  the  p-side  into  four  regions  similar  to  the  double 
injection  with  trapping  case. 


V(volts) 

Figure  22.  Forward  I-V  characteristics  for  alloyed  germanium  diodes  (Ref.  14). 
The  graph- displays  the  JaV^  initial  characteristic  expected 
before  second  breakdown. 
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the  applied  polarities).  Solving  one  side  of  this  problem  for  the  general 
case  provides  a  solution  for  the  other  side. 


Page 34  provides  the  full  problem  description  and  general  solution  for 
this  problem.-  The  threshold  voltage  for  breakdown  is 


Vtw  =  ■  ■  =  - 


«  lOV 


2/Ap  ifipTp 

A  smlllar  expression  holds  for  the  N-side 


(204) 


Vth  =  7V 


ilinTn 


(205) 


where  L  is  the  length  of  each  bulk  side,  and  adds  for  the  devices  total 
potential  difference. 


The  minimum  voltage  reached  after  breakdown  is 

?(<»)  VMp^p/  \Mn»-n/S(a) 


(206) 


with  the  associated  minimum  current  density 


Jm  =  h{a) 


^SrL 


(207) 


As  demonstrated  by  Equations  204  and  205,  a  situation  arises  where  the 
transit  time  equals  the  lifetime  of  the  carriers.  Also  known  as 


transit  time  »  Lifetime  = 


12 


VthI*p 


(208) 


a  transition  from  single  carrier  injection  to  double  carrier  injection.  The 
ratio  of  these  quantities  yields  the  concept  of  gain. 


“^dbllnj  _  t  _  ^  _ 


7. 


Slinj 


id  t 


gain 


(209) 


When  t  is  equal  to  or  less  than  t  the  gain  and  double  injection 
dominates. 
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To  solve  for  Vth*  and  for  this  diode  we  have  a  =  4  on  the  p-side 
and  a  *  4  on  the  n-side. 

Vth  (over  total  length)  =  17  V  . 

Vm  (over  total  length)  =  8.63  V 

Jni  (over  total  length)  =  26.7  A/cm^  *»  0.2  A 

If  thermal  effects  were  allowed  into  the  model,  different  results  would 
certainly  occur.  Starting  the  simulation  at  room  temperature  would  certainly 
ionize  all  impurities,  while  the  electrical  heating  of  the  device  would 
quickly  make  the  material  intrinsic  and  lower  carrier  mobility  and  decrease 
lifetime.  This  would  decrease  the  value  of  Vth»  and  Jp,.  A  number  of 
solutions  are  then  possible  dependent  on  the  rate  of  heating  and  heat  loss. 

Use  of  Ridley's  equations  would  yield  the  filament  cross  section  near  60 
percent  of  the  device's  area  upon  reaching  Vm  and  jm  ignoring  thermal  effects. 
Adding  thermal  effects  would  no  doubt  increase  the  filament  cross  section. 

Determining  how  the  simulation  would  progress  after  filament  formation 
with  heating  would  quickly  enter  the  realm  of  solid  state  plasma  physics  which 
is  beyond  the  scope  of  this  effort.  The  only  definitive  statement  that  can  be 
made  is  sustaining  operation  after  filament  formation  wi.th  heating  will  lead 
to  device  degradation  via  thermal  shock. 

PN  JUNCTION— REVERSE  BIAS 

Under  reverse  bias,  the  diode  considered  here  behaves  quite  differently. 
This  analysis  will  concentrate  on  the  aspects  that  occur  after  avalanche  has 
begun . 

The  depletion  region  continues  to  grow  as  the  applied  potential  increases. 
Thus,  close  attention  will  not  have  to  be  paid  to  the  bulk  regions.  The 
depletion  region  has  two  distinct  areas  with  different  characteristics. 

The  avalanche  region  is  basically  neutral  from  the  electron-hole  pairs 
being  created  from  impact  ionization.  The  region  involved  in  avalanche  is 
less  than  one-tenth  the  length  of  the  depletion  region  for  a  two-sided  abrupt 
junction  (Ref.  15)  (Fig.  23).  The  remaining  portion  of  the  depletion  region 
consists  of  traps  that  are  partially  filled.  As  carriers  from  the  avalanche 
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AREA  UNDER  AVALANCHE 


Figure  23.  Depiction  of  the  depletion  region  of  a  reverse  biased  abrupt  PN 
junction.  The  shaded  area  is  under  avalanche  injecting  carriers 
into  the  drift  regions  filled  with  traps  (not  to  scale). 

region  enter  the  recombination  zone  space  charge  builds  up  from  the  inability 
of  the  majority  of  carriers  to  transit  as  region  without  being  trapped. 

The  injection  from  the  bulk  regions  is  very  small  and  thermal  excitation 
is  not  adequate  to  relax  the  space  charge.  In  addition,  as  the  applied  voltage 
increases,  the  depletion  region  enlarges  with  the  square  root  of  the  voltage. 

The  problem  can  be  divided  up  as  noted  in  Fig.  24.  The  region  from  x  =  0 
to  Xi  is  the  avalanche  region  and  is  only  a  fraction  of  the  width  x.  As  the 
applied  voltage  increases  the  fraction  of  x  used  by  this  region  is  still  quite 
small  until  very  high  fields  (over  10®  V/cm)  are  formed.  The  region  from  xj 
to  X2  is  a  transition  region  that  varies  in  size  dependent  on  the  amount  of 
thermally  excited  carriers  and  the  avalanche  regions  size.  From  X2  to  Xp  the 
trapped  positive  space  charge  dominates  current  flow.  Basically,  a  recombina¬ 
tion  barrier  exists  that  current  flow  (from  x  =  0  to  Xp)  must  overcome.  Single¬ 
carrier  injection  is  the  prime  mode  of  current  flow. 

Application  of  the  Regional  Approximation  Method  to  this  type  of  problem 
was  accomplished  in  solving  the  trap-filled  limit  problem  (p.  9).  The  neutral 
region  expands  with  increasing  voltage  until  the  trapped  space  charge  is  over¬ 
come,  leading  to  double-carrier  injection  and  CCNR. 
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figure  24.  N  side  of  the  depletion  region  (not  to 
scale). 


All  three  regions  expand  with  applied  voltage,  yet  the  neutral  region 
will  not  expand  to  fill  the  space  charge  region  until  extremely  high  fields 
are  attained".  This  is  only  probable  for  devices  at  liquid  nitrogen  tempera¬ 
tures  where  effects  due  to  intrinsic  carrier  increases  are  minimal.  At  or 
near  room  temperature  the  increasing  quantity  of  thermal  carriers  would  be  the 
primary  cause  of  space  charge  relaxation  leading  to  double  injection  (i.e., 
the  point  where  the  gain  is  greater  than  or  equal  to  one). 

'Starting  with  the  operating  trap-filled-limit  regime,  a  method  of  attach¬ 
ment  must  be  developed  which  allows  the  use  of  changing  lengths  dimension  and 
heating.  An  examination  of  p.  14  shows  the  critical  current  and  voltage  are 


2e^NtnL  „  4«Ar,l2 

"  ''"-Sir 


(210) 


for  a  crystal  of  length  L.  The  depletion  layer  width,  x^,  is  defined  in  terms 
of  the  applied  voltage  as 


(kiVii[ND  +  N,)\^ 

\  ) 


(211) 


where 

Np  +  N^  _L 
NpNa  "  Nt 

The  equations  noted  in  Equation  210  apply. to  the  insulator  part  at  the  depletion 
region  after  avalanche  occurs.  This  region  occupies  around. 90  percent  of  the 
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depletion  region.  The  voltage  across  will  be  the  voltage  across  the  ava¬ 
lanche  region  and  the  voltage  across  the  insulator. 

The  problem  is  simplified  initially  by  ignoring  thermal  effects  and  hold¬ 
ing  the  avalanche  region  at  a  constant  length  of  0.1  Xp,  so  that 

=  (212) 

where 

Em  =  2.25  X. lO^  V/cm 
Xp  =  1.4  X  10“^  cm 
Vr  =  337  V 

From  which  we  find  i|»aval  ’  16.65  v  and  the  remaining  potential  drop  across  the 
Insulator  part  of  the  depletion  region  is  about  152  V.  Solving  for  the  voltage 
necessary  from  Equation  197,  V^r  is  about  362  V.  To  decrease  the  value  of 
L  must  be  small  compared  to  the  avalanche  region  to  overcome  Thus, 

Xg  »  xp  -  Xa  before  is  achieved  (for  this  particular  example  Xg  'v  .9  Xp 
before  V^r  is  overcome,  thus  thermal  effects,  are  rejoined  to  aicl  the  device). 

If  thermal  effects  are  allowed  in  an  isothermal  situation,  it  is  apparent 
that  increasing  intrinsic  carrier  concentration  will  lower  in  Equation  210 
by  filling  these  traps.  Thus,  V^-t  should  become. lower  with  higher  temperature. 
The  rate  of  heat  absorbed  and  lost  will  determine  when  is  lowered  for  a 
particular  Vr. 

This  thermal  assist  mechanism  can  have  a  lesser  effect  if  reemission  of 
carriers  is  allowed  in  the  insulator.  Also,  the  addition  of  other  impurities 
can  affect  lifetime  and  decrease  or  increase  the  chance  for  a  thermal  assist. 

Solving  these  equations  with  all  effects  would  require  a  great  deal  of 
knowledge  about  the  semiconductor.  Pertinent  data  left  out  of  the  equations 
can  result  in  a  large  error.  Effects  due  to  very  high  fields  such  as  changing 
mobilities,  have  not  been  considered  and  are  beyond  the  scope  of  this  report. 

By  examining  the  PN  diode  in  forward  and  reverse  bias  the  Regional 
Approximation  method  allows  a  definite  insight  to  determine  what  precedes 
CCNR. 
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Both  situations  examined  for  a  diode  result  in  relaxation  of  an  accumu¬ 
lated  space  charge  by  increased  single  carrier  injection.  As  single-carrier 
injection  dissipates  the  space  charge,  double-carrier  injection  follows  with 
CCNR  and  filament  formation. 

The  particular  examples  studied  here  show  thermal  effects  aid  the  elec¬ 
tronic  process  by  relaxing  space  charge  build-up.  Thus,  the  time  history  of 
use  and  variables  that  effect  the  thermal  characteristics  of  a  device  do  play 
a  role  in  breakdown  for  these  particular  examples.  Other  device  differences 
(such  as  graded  or  exponential  junctions)  will  require  study  on  a  case-by-case 
basis  to  discern  the  transport  mechanism  but  will  still  follow  the  same  general 
trends  as  these  examples.  . 

It  is  important  to  note  that  thermal  effects  alone  do  not  cause  or  initiate 
the  occurrence  of  CCNR.  Thermal  effects  assist  the  initial  conditions  or  relax 
space  charge,  but  it  is  the  electronic  transport,  or  transport  gain  that  causes 
CCNR  and  its  accompanying  effects. 

PIN  DIODE- FORWARD  BIAS 

This  type  of  diode  is  used  for  high  voltage  applications  primarily.  The 
extra  I  diffusion  provides  an  additional  expansion  region  for  the  depletion 
region  under  high  reverse  and  forward  bias.  The  diode  considered  is’  detailed 
below  and  in  Fig.  25. 

P=N=lxl0^5  cm" 3,  abrupt  diffusions 

I  (n-type)  -lx  10^^  (-^-3 

cylindrical  geometry 

radius  *  1  x  lO"^  cm 

Tn  *  Tp  *  1  X  10-6  s 

un  =  1350  cm2/V-s 

Wp  *  480  cm^/V-s 

Built-in  potential  (PI  junction)  =  0.397  V 
Unbiased  depletion  layer  width  (PI  junction)  =  7.19  x  10"“+  cm 
Built-in  potential  (IN  junction)  =  0.119  V 
Unbiased  depletion  layer  width  =  3.94  x  10’**  cm 
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Figure  25. 


PIN  diode  under  consideration,  not  to 
scale. 


Under  high  forward  bias  the  depletion  regions  of  a  PIN  diode  will  be 
reduced  in  size  in  a  manner  similar  to  a  PN  diode.  The  IN  junction  will  be 
negated  first  by  the  applied  voltage.  The  PI  junction  will  eventually  be 
negated,  but  not  till  a  potential  develops  across  the  P  and  I  regions.  This 
potential  gradient  creates  an  electric  field  and  the  initiation  of  carrier 
drift.  Holes  are  injected  from  the  P  region  into  the  I  region,  while  electrons 
are  injected  from  the  IN  Junction.  Region  I  is  characterized  as  an  insulator 
with  fully  ionized  donor  centers,  which  represent  a  recombination  barrier  to 
electrons. 


Initially,  the  J-V  characteristics  will  reflect  a  space  charge  limited 
current  where  J  «  (Ref.  2).  As  the  current  continues  to  increase,  and 
direct  recombination  of  carriers  occurs  the  current  transport  characteristics 
transform  to  a  semiconductor  injected  plasma.  The  J-V  characteristics  then 
conform  to  J  «  V^  (Ref.  2)  as  noted  in  Fig.  13.  As  before,  the  threshold 
voltage  is  the  same  on  a  forward  biased  PN  diode. 


Vth  = 


2^^nrn 


where  the  average  transit  time  equals  the  lifetime  of  electrons  in  the  I 
region  (the  effects  on  forward  voltage  are  seen  in  Ref.  16.  For  the  PIN  diode 
described,  the  across  the  I  region  is  about  9.2  mV.  The  voltage  will  drop 
at  this  point  to  about  4  mV  across  Region  I  while  the  current  continues  to 
increase.  Heating  rapidly  destroys  the  device  after  CGNR  occurs. 
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These  effects  assume  that  the  lifetime  and  mobility  (i.e.,  scattering 
mechanisms)  are  not  altered  by  injection  level,  and  space  charge  build-up  is 
confined  to  very  near  the. junctions.  Under  realistic  conditions  the  relaxa¬ 
tion  of  these  requirements  will  increase 

PIN  DIODE— REVERSE  BIAS 

Under  reverse  bias,  the  PIN  diode  resembles  a  long  diode.  At  low  and 
moderate  voltages  the  diode  exhibits  a  fairly  low  current  J  «  V  tending  toward 
J  «  V^.  The  depletion  region  exists  primarily  in  Region  I,  i.e.,  the  insula¬ 
tor  region.  As  the  diode  nears  avalanche,  the  depletion  region  normally  fills 
the  insulator  entirely.  As  avalanche  occurs  the  entire  insulator  becomes 
filled  with  electron-hole  pairs  near  the  saturated  drift  velocity.  Little  of 
the  depletion  region  extends  into  the  P  region  and  practically  nothing  extends 
into  the  N  region. 

If  the  avalanche  coefficients  are  allowed  to  be  equal,  the  avalanche 
breakdown  voltage  takes  the  form 


Wb 


i9l8  V 


(213) 


where  _ 

b  =  1.65  X  106  v/cm 
a  =  1.6  X  10®  V/cm 
w  =  insultator's  width  in  cm 

With  an  E„,„  =  1.83  x  10®  V/cm  constant  across  the  insulator.  The  electron- 
hole  pairs  produced  in  the  insulator  separate,  holes  leave  the  PI  interface 
via  diffusion  and  electrons  leave  the  IN  interface  by  diffusion.  A  small 
portion  of  the  depletion  region  extends  into  the  P  region,  but  none  penetrates 
the  N  region.  Therefore,  very  few  carriers  trespass  the  junctions  in  large 
numbers  since  diffusion,  not  drift,  dictates  the  current  level. 

The  depletion  region  extends  about  5  x  10"®  cm  into  the  P  region  with  a 
voltage  drop  of  about  4.5  V.  The  only  way  to  open  the  N  junction  for  a  large 
carrier  flow  is  by  increasing  the  potential  till  the  N  Junction  reaches  ava- 
lance  allowing  the  space-charge  accumulating  in  the  insulator  to  relax  (Fig. 
26). 
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Figure  26.  Change  of  the  electric  field  with 
position  as  the  potential  across 
the  device  increases.  The  field 
mimics  the  charge  distribution 
(not  drawn  to  scale). 

The  distance  the  depletion  ragion  must  penetrate  into  the  P  and  N  regions 
must  reach  at  least  one  debye  length  if  avalanche  is  to  be  sustained  (for  the 
avalanche  multiplication  to  approach  infinity). 

The  extrinsic  debye  length,  Lq  is 

/  ''  =  ’'7  -5 

'  ^  '  '  qlNp  -  n  '  '  ’  5  ■<  1°  cm 
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To  attain  this  distance  on  the  N-side  requires  a  potential  drop  of  about  45  V. 
Thus  the  total  voltage  is  50  V  (PI  Junction)  plus  918  V  (for  the  insulator) 
and  45  V  (at  the  IN  junction)  a  total  of  1013  V.  At  or  near  this  voltage  the 
N  region  junction  will  begin  to  avalanche  injecting  large  numbers  of  carriers 
into  the  insulator  and  bulk  regions  passing  a  large  current  and  reducing  the 
insulators  space-charge  and  field  strength.  The  diode  has  transferred  from 
single  carrier  injection  dominated  current  flow  (SCL  current  flow)  to  double¬ 
carrier  injection  in  a  trap-free  insulator— also  called  the  semiconductor 
j  oe  regime.  Fi lamentation  follows  the  CCNR  quite  quickly.  Experimental 
proof  verifying  the  preceding  description  is  given  in  Refs  15  and  17- 

If  thermal  affects  are  allowed  to  enter  into  the  situation  several  changes 
occur.  Avalanche  coefficients  will  decrease  in  magnitude  due  the  phonon  inter¬ 
action  with  free  carriers.  The  intrinsic  carrier  concentration  will  increase 
rapidly  (doubling  for  every  11  rise  above  room  temperature)  to  effect  tne 
occupation  of  donor  sites  in  the  insulator.  As  the  sites  fill,  the  material 
becomes  intrinsic  decreasing  the  mobility  and  possibly  changing  the  carrier 
lifetime.  The  increased  carrier  density  in  the  bulk  regions  will  encourage 
farther  diffusion  and  heighten  the  overall  current  passed  by  the  diode.  To 
achieve  avalanche  in  the  bulk  regions  near  the  junctions  will  therefore  require 
a  higher  potential.  The  voltage  needed  to  breakdown  the  diode  will  therefore 
be  higher. 

The  PIN  diode  shows  many  of  the  same  conditions  that  precedes  second 
breakdown  in  the  PN  diode.  Under  reverse  bias,  the  initiation  of  conditions 
for  second  breakdown  require  avalanche  breakdown  in  the  bulk  regions  of  the 
junctions.  This  is  followed  by  space  charge  relaxation  in  the  insulator  and 
double  injection. 

BIPOLAR  JUNCTION  TRANSISTORS— SINGLE  DIFFUSED 

The  Bipolar  Junction  Transistor  (BJT)  is  the  next  step  in  the  evaluation 
of  diode  technology.  Today’s  BJTs  are  fabricated  with  many  different  types  of 
diffusion  processes  and  numerous  architectures.  Application  of  the  Regional 
Approximation  method  to  the  transistor  has  been  performed  (Ref.  IS)  and  wir 
not  be  repeated  here.  Instead,  the  single  diffused  somewhat  idealized  transis¬ 
tor  will  be  analyzed  for  special  situations  where  large  voltages  are  applied 
in  various  ways  to  achieve  unstable  device  operation.  This  can  then  be 
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examined  as  special  cases  of  unstable  current  states  (CCNR)  as  shown  for  the 
PN  and  PIN  diodes.  Cases  Involving  double-diffused  BJT  and  field  effect 
transistors  will  not  be  examined.  The  reader  will  find  after  some  considera¬ 
tion  that  these  devices  are  also  special  cases  of  the  PN  or  PIN  diodes. 

The  generic  BJT  examined  here  will  not  be  quantified  numerically  with 
general  mathematical  criteria  for  unstable  operation.  The  numerous  types  of 
BJTs  and  prolific  architectures  that  are  available  often  times  dictate  how  the 
transistor  parameters  (a,  6,  spreading  resistance  etc.)  are  related  at  high 
bias.  Therefore,  general  cases  would  have  to  be  developed  for  each  particular 
BJT  type  and  architecture,  which  Is  beyond  the  scope  of  this  report. 

The  device  to  be  examined  Is  a  P-N-P  single  diffused  BJT  as  shown  in  Fig. 
27.  Then,  the  BJT's  most  important  dimension  Is  base  width  and  contact  peri¬ 
phery.  To  achieve  transistor  action,  the  amplification  of  a  signal,  the  emitter 
must  be  close  enough  to  the  collector  to  Influence  it.  If  the  base  is  wide 
compared  to  the  carrier  diffusion  length,  excess  carriers  injected  from  the 
emitter  cannot  reach  the  collector  and  Influence  It.  Therefore,  the  base 
width  Is  on  the  order  of  a  diffusion  length.  The  base  width  directly  influences 
a  factor  called,  the  base  transport  factor,  s.  Beta  represents  the  ratio  of 
minority  carrier  current  reaching  the  collector  to  the  actual  minority  carrier 
current  emitted.  As  such.  It  Is  also  a  measure  of  the  carrier  transit  time 
divided  by  the  lifetime.  The  value  for  s  is  usually  slightly  less  than  uni  ty. 

Another  factor  used  to  qualify  the  BJT  is  6,  the  collector  efficiency. 

This  represents  the  ratio  of  collected  to  Incident  current  on  the  collector 
and  is  normally  less  than  unity.  The  final  value  of  this  variable  is  dependent 
on  the  value  of  the  surface  recombination  rate. 

The  symbol  y  is  called  the  emitter  injection  efficiency  and  represents  the 
ratio  of  minority  current  to  total  emitter  current.  The  value  of  y  is  normally 
less  than  unity. 

All  three  factors  listed  above  can  be  mtl tipi  led  together  to  yield  the 
factor  og  known  as  the  current  ampi if 'cation  factor.  It  is  aouroxlmately  equal 
to  the  negative  of  the  total  collector  current  minus  the  collector  leakage 
current  divided  by  the  total  emitter  current. 

a, (21^5) 
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Figure  27.  Schematic  diagram  of  a  pnp  transistor  under  no  external  bias  (not 
drawn  to  scale).  All  junctions  are  abrupt  diffusions. 

Note  that  under  any  bias  condition,  either  (I^  -  or  will  be  negative 
by  definition  making  og  positive.  Many  aspects  of  transistor  operation  can  be 
explained  in  terms  of  these  variables  including  aspects  of  electronic  instabi¬ 
lity  in  the  four  cases  to  be  considered.  '  * 

COLLECTOR  BREAKDOWN 

A  p-n-p  transistor  is  illustrated  in  Fig.  28  in  the  common  base  configura¬ 
tion 

The  first  extreme  of  operation  to  be  considered  is  avalanche  breakdown  of 
the  collector  Junction  with  the  emitter-base  loop  open.  This  is  nothing  more 
than  simple  diode  breakdown  already  outlined  in  an  earlier  section. 

In  terms  of  the  transistor  variables  outlined  before,  the  actual  number 
of  minority  carriers  collected  becomes  much  larger  than  the  number  of  minority 
carriers  that  actually  strikes  the  collector.  This  is,  of  course,  due  to 
avalanche  breakdown  of  the  collector. 

Numerically,  a  is  determined  by  the  multiplication  coefficient  of  avalanche. 
This  has  been  defined  experimentally  as 


n  ■  2  fo  4  (216) 
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(Figure  28.  A  pnp  transistor  in  the  common  base  configuration.  The  emitter  is 

forward  biased  with  respect  to  the  base.  The  collector  is  reverse 
biased  with  respect  to  the  base.  Drawing  is  not  to  scale. 

1  As  5  tends  toward  infinity,  so  does  ag  in  Equation  215.  An  electrical  instabil- 

j  ity  will  occur  as. described  previously  for  the  PN  diode  in  reverse  bias.  But, 

I  two  limitations  apply  specifically  to  transistors  in  this  scenario.  The  col¬ 

lector  depletion  region  cannot  expand  beyond  the  limit  of  the  base  region 
without  support  of  an  external  potential.  Since  the  base  electrode  is  normally 
located  around  the  periphery  of  the  base  diffusion,  the  depletion  region 
expands  preferentially  toward  those  regions.  If  the  depletion  region  touches 
the  electrodes  the  space  charge  resistance  (due  to  SCL  current  flow)  will  be 
dissipated  initiating  double  carrier  injection.  The  inclusion  of  surface 
recombination  or  spreading  resistance  makes  punch-through  to  an  electrode  much 
harder.  Heating  can  bring  the  material  into  the  intrinsic  regime  relaxing  the 
accumulated  space  charge  and  possibly  initiating  solid  state  plasma  formation. 

EMIHER— COLLECTOR  BREAKDOWN 

This  situation  is  depicted  in  Ffg.  29  with  the  emitter  junction  forward 
biased  and  the  collector  junction  reverse  biased.  Initially,  the  base  i5  to 
be  open. 

Since  the  emitter  is  forward  biased,  there  is  a  minority  current  into  the 
base.  If  Vq£  is  greater  then  a  few  tenths  of  a  volt,  the  emitter  is  Injecting 
a  current  into  the  base  that  is  greater  than  the  collector's  reverse  saturation 
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Figure  29.  Set-up  to  study  collector-emitter  breakdown  (not  drawn  to 
scale) . 


current.  Thus,  the  current  through  the  collector  is  larger  than  This 

current  is  plus  Ij,g  or 


as  r  1 

as  =  - - « 

1  —  ag  Id 


icEo  =  Icbu(1-^e)  (217) 

As  Vgg  increases,  the  width  of  the  collector  depletion  region  expands  into 
the  base  increasing  g  and  5.  also  increases  and  conversely  y.  These 
increases  boost  05  and  Ig  until  avalanche  is  initiated.  The  increased  quantity 
of  ca-rriers  in  the  collector  junction  means  avalanche  is  initiated  at  a  lower 
field  strength  than  BV^gg. 

Further  increases  in  Vgg  result  in  the  collector  depletion  region  reaching 
the  emitter  junction  which  provides  carriers  to  relax  the  accumulated  space 
charge  initiating  double  injection  and  CCNR.  This  situation  is  quite  similar 
to  breakdown  in  a  PIN  diode. 


In  terms  of  the  transistor  variables,  5  does  not  increase  to  infinity 
with  higher  Vgg  until  it  is  closer  to  avalanche.  But  B  and  y  increase  to  high 
values  since 


Ofl  ,  1 

- m  I.'  ~  *  — - - 

1  -  oq  10  I  - 


(218) 
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The  current  Iq,  rises  dranatlcally.  As  punch- through  occurs,  the  space  charge 
relaxes,  initiating  double  injection  and  CCNR  (a^  will  increase  very  rapidly  as 
Og  increases  to  infinity  with  punch- through)  as  shown  in  the  plot  of 
Fig.  30. 

Now,  allow  the  base  to  be  connected  back  to  the  circuit  in  Fig.  29.  The 
value  Ig  is  not  zero  and  is  able  to  drain  carriers  out  of  the  base.  The 
inmediate  consequence  of  Ig^^O  is  that  decreases  from  the  last  case 
dependent  on  the  value  of  og  at  low  level  Injection. 

As  Vq£  is  increased  the  last  scenario  is  repeated  with  the  following  added 
complications.  Since  IgfO.  current  is  flowing  toward  the  periphery  of  the 
base  to  the  electrical  contacts.  There  is  little  resistance  to  carrier  flow 
perpendicular  to  the  junctions,  but  parallel  flow  meets  considerable  resistance 
because  of  the  vector  quality  of  carrier  diffusion  and  drift  in  the  base  (Ref. 
19).  The  exact  amount  of  resistance  is  determined  by  the  base  width  and  overall 
geometry  (Ref.  20).  The  net  result  of  this  spreading  resistance  is  to  make 
the  center  of  the  emitter  less  forward  biased  then  the  edges.  Thus,  current 
flow  is  higher  at  the  edges  than  the  center  and  the  Early  Effect  serves  to 
further  enhance  the  overall  potential  build-up.  The  value  of  Vc  increase 
past  BV^q  until  the  potential  across  the  accumulated  space  charge  and  spreading 
resistance  is  large  enough  to  allow  the  carrier  transit  time  to  equal  the  car¬ 
rier  lifetime. 

The  resulting  Iq  versus  V^.  curve  will  behave  as  shown  in  Fig.  31.  Curve 
B  in  Fig.  31  shows  a  CCNR  aspect  that  results  even  if  avalanche  has  already 
occurred  at  the  collector  junction.  Such  activity  is  specific  to  the  BJT  and 
is  dependent  on  geometry  and  the  external  circuit.  This  type  of  activity  has 
been  demonstrated  in  Reference  21. 

ELECTRICAL  INSTABILITY  DURING  COMMON  EMIHER  MODE  CUTOFF 

The  situations  described  previously  cover  most  extreme  forms  of  operations 
that  can  be  discerned  for  a  BJT  used  as  a  switch  or  amplifier.  One  case  does 
remain  where  both  emitter  and  collector  are  reverse  biased  simultaneously. 

This  occurs  when  a  BJT  is  used  in  the  common  emitter  set-up  shown  in  Fig.  32. 

As  Vq£  and  V^g  increase,  the  base  becomes  modulated  by  the  expanding 
depletion  regions.  As  the  voltage  reaches  avalanche  for  both  junctions  the 
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Figure  30.  General  I-V  curves  showing  breakdown  for 
the  collector-base  compared  to  the  col¬ 
lector  emitter  breakdown  (not  to  scale). 


Figure  31.  General  I-V  curves  for  coll ector-emi tier 
breakdown.  Curve  A  is  collector- 

emitter  avalanche  breakdown  with  no 
resistance  in  the  base.  ^ 

possible  due  to  space  charge  build-up, 
base  spreading  resistanre.  and  some 
external  resistive  oa''i  net  ‘o  scale 
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Figure  32.  A  pnp  transistor  in  the  common  emitter  setup  (not  to  scale).  Both 
junctions  will  be  reverse  biased  to  create  a  cutoff  situation. 


current  injected  into  the  base  will  be  affected  by  base  spreading  resistance. 
This  will  reduce  the  potential  near  the  center  of  each  junction  forcing  cur¬ 
rent  flow  to  the  edges.  The  magnitude  of  this  effect  will  be  greatly  dependent 
on  geometry  and  the  amount  of  current  injected  into  the  base.  If  the  depletion 
regions  expand  into  the  base  and  merge  we  have  the  breakdown  problem  of  the" 

PIN  diode.  Space  charge  is  relaxed  and  double  carrier  injection  results.  The 
Ij.  -  curve  will  behave  as  shown  in  Fig.  33. 

For  curve  A,  the  base  offers  no  resistance  to  current  flow  resulting  in 
avalanche,  punch-through,  and  double  carrier  injection.  In  terms  of  the 
transistor  variables,  becomes  very  large  and  tends  to  infinity  once  punch- 
through  occurs. 


For  curve  B,  spreading  resistance  due  to  the  removal  of  carriers  from  the 


base  results  in  an  increase  in  before  avalanche  can  occur.  This  curve  will 


also  depend  on  the  device  geometry  and  external  circuit.  This  effect  may  have 
been  observed  in  Reference  21  numerous  times.  Also,  curvature  of  the  emitter- 
base  deplet  on  region  (due  to  the  emitter  area  being  smaller  than  the 
collector  area)  results  in  preferential  avalanche  near  the  regions  of  highest 
curvature  while  the  collector-base  junction  will  avalanche  homogeneously.  The 
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Figure  33.  General  I-V  curves  for  breakdown 
with  the  transistor  in  cutoff. 
Curve  A  is  for  no  resistance  in 
the  base  loop.  Curve  B  is  for 
resistance  in  the  base  loop.  The 
amount  of  resistance  will  deter¬ 
mine  the  shape  of  curve  B  (not 
drawn  to  scale). 


emitter-base  avalanche  is  possibly  current  limited  by  surface  recombination 
and  proximity  to  the  base  contacts.  This  will  not  allow  s  or  y  to  increase 
with  V^g,  and  may  very  well  decrease  with  increasing  V^g.  It  is  even  possible 
that  surface  breakdown  will  limit  V^g  before  homogeneous  avalanche  can  occur. 


If  the  limitations  of  V^g  are  ignored  for  now,  the  increased  values 
needed  before  breakdown  can  occur  are  explainable  by  resistance  effects. 
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CONCLUSIONS 

The  use  of  the  Regional  Approximation  method  allows  in-depth  analysis  of 
devices  under  high  level • injection.  Basic  scenarios  explaining  how  electron¬ 
ically  unstable  phenomena  can  occur  have  been  detailed  and  applied  to  a  simple 
PN  and  PIN  diode.  Application  to  transistors  were  detailed  for  cases  of 
extreme  operation. 

An  exhaustive  application  was  not  performed  for  all  types  of  device  dif¬ 
fusions,  materials,  or  architectures.  Instead,  basic  electronic  transport 
problems  in  insulators  and  semiconductors  was  developed  to  show  how  injection 
level  plays  a  significant  role  in  determining  J-V  curves.  These  cases  were 
shown  to  apply  to  device  physics  for  a  few  major  scenarios. 

A  vigorous  connection  was  developed  between  the  transition  of  single-  to 
double-carrier  injection  and  the  initiation  of  current  controlled  negative 
resistance.  This  provides  a  basis  for  understanding  the  formation  of  current 
filamentation  as  a  result  of  an  electronically  unstable  situation. 

Other  methods  of  obtaining  electronically  unstable  states  relates  to  the 
quantum  mechanical  aspects  of  the  material  itself  as  elucidated  in  the  com¬ 
panion  report,  "Analysis  of  Solid  State  Plasma  Formation  in  Semiconductor 
Components,"  AFWL-TR-85-116. 
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